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1.  INTRODUCTION 


1.1  GOALS  AND  OBJECTIVES 

The  objective  of  this  dissertation  is  to  describe  wave  propagation  in  a  fluid-loaded, 
homogeneous,  transversely  isotropic,  elastic  cylinder  of  arbitrary  thickness  and  to  in¬ 
vestigate  the  effects  of  fluid  loading  on  the  cut-off  frequencies  and  phase  velocities  of 
various  axial  shear  and  plane  strain  modes. 

In  particular,  the  goal  is  to  investigate  the  cut-off  frequencies  and  phase  velocities 
of  the  flexural  modes  (n=l,  where  n  is  the  circumferential  wavenumber)  of  wave  prop¬ 
agation  in  the  coupled  system  consisting  of  afiber-reinforced  cylinder  that  is  filled  with 
and  immersed  in  an  inviscid  fluid.  To  assist  in  the  interpretation  of  the  numerical  re¬ 
sults,  an  analysis  of  the  n=l  modes  of  wave  propagation  in  fluid-loaded  isotropic  cyl¬ 
inders,  with  properties  of  steel  and  of  soft  (linear)  rubber,  is  also  conducted. 

Since  the  primary  interest  is  in  the  propagating  modes,  branches  extending  into  the 
imaginary  and  complex  wavenumber  domains  are  not  included  in  this  study.  However, 
the  attenuation  of  these  modes  due  to  the  external  fluid  loading  can  be  obtained  directly 
from  the  imaginary  part  of  the  complex  frequency. 

1.2  BACKGROUND 

1.2.1  Basic  concepts  in  wave  propagation 

The  following  definitions  are  provided  as  background  information.  The  text  by 
Achenbach  (1990)  provides  a  comprehensive  study  on  wave  propagation  in  elastic  sol¬ 
ids. 


1 


1.2. 1.1  Phase  velocity 


A  traveling,  mechanical  wave  in  one  dimension  is  defined  by  an  expression  of  the 
type  f  =  f  (cot  -  kz) ,  where  f,  as  a  function  of  the  spatial  coordinate  z,  time  t,  wave- 
number  k,  and  radian  frequency  co ,  generally  denotes  a  displacement,  a  particle  velcx;- 
ity,  or  a  stress  component.  The  argument  ©t  -  kz  is  the  phase  of  the  wave  function. 
Points  of  constant  phase  are  propagated  with  the  phase  velocity 


c  =  ©/k.  (1.1) 

In  an  inviscid,  compressible  fluid,  the  phase  velocity  is  the  velocity  of  sound  in  the  fluid 
(or  the  acoustic  phase  velcxiity). 

If  c  =  ©/k  is  complex,  the  wave  function  is  attenuated.  The  phase  velocity  is  the 
real  part  of  c.  The  attenuation  is  defined  by  the  imaginary  part  of  c.  Normally,  the  at¬ 
tenuation  is  written  as  an  exponential  decay  as  a  function  of  time  t  or  of  the  spatial  co¬ 
ordinate  z. 

1.2. 1.2  Dispersion 

If  the  phase  velocity  depends  on  the  wavenumber,  the  system  is  said  to  be  disper¬ 
sive.  Dispersion  cx:curs  in  inelastic  (dissipative)  materials  and  in  elastic  waveguides, 
where  a  waveguide  is  any  extended  body  with  a  cross  section  of  finite  dimensions  (such 
as  an  infinite  plate  and  an  infinitely  long  cylinder). 

1.2. 1.3  Group  velocity 

Group  velocity  is  the  velocity  of  energy  transmission  in  the  system.  Group  velocity 
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is  defined  as 


=  dco/dk.  (1.2) 

e 

Using  CO  =  ck ,  we  can  write  the  group  velocity  as 

c  =  c  +  k  dc/dk.  (1.3) 

e 

If  there  is  dispersion,  the  group  velocity  is  different  from  the  phase  velocity, 

1.2. 1.4  Frequency  equation 

The  frequency  equation  is  an  expression  relating  co  to  in  waveguides.  The  fre¬ 
quency  equation  is  generally  derived  from  the  equations  of  motion  and  the  boundary 
conditions  of  the  system. 

1.2. 1.5  Branches  and  dispersion  curves 

Solution  of  the  frequency  equation  yields  an  infinite  number  of  continuous  curves 
called  branches.  A  branch  displays  the  relationship  between  wavenumber  and  frequen¬ 
cy  for  a  particular  mode  of  propagation.  If  the  branch  is  dispersive,  it  is  sometimes 
called  a  dispersion  curve. 

1.2. 1.6  Frequency  spectrum 

A  collection  of  branches  is  the  frequency  spectrum  of  the  system.  A  complete  fre¬ 
quency  spectrum  includes  branches  extending  into  the  imaginary  and  complex  wave- 
number  domains. 
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1.2. 1.7  Cut-off  frequency 


Cut-off  frequencies  are  the  frequencies  corresponding  to  a  zero  wavenumber.  For  a 
particular  mode,  the  cut-off  frequency  is  the  frequency  where  the  wavenumber  of  the 
branch  changes  from  real  to  imaginary  (or  complex). 

1.2.2  Transversely  isotropic  materials 

Many  natural  and  artificial  materials  are  transversely  isotropic.  For  example,  hex¬ 
agonal  crystals,  including  beryllium,  cadmium,  magnesium,  titanium  and  zinc,  are 
transversely  isotropic.  Ice  is  also  transversely  isotropic. 

A  unidirectional,  fiber-reinforced  composite  material  is  transversely  isotropic.  This 
type  of  material  is  highly  anisotropic  because  the  stiffness  and  strength  in  the  fiber  di¬ 
rection  are  of  the  order  of  the  values  for  the  fiber,  and  are  thus  very  high,  and  the  stiff¬ 
ness  and  strength  transverse  to  the  fiber  direction  are  of  the  order  of  the  values  for  the 
matrix,  and  are  thus  much  lower.  Carbon  and  graphite,  which  are  often  used  in  fiber- 
reinforced  materials,  are  themselves  transversely  isotropic. 

1.2.2. 1  Constitutive  behavior  of  transversely  isotropic  materials 

In  its  most  general  form  (here,  and  in  the  following,  we  use  the  Cartesian  coordinate 
system),  Hooke’s  law  can  be  written  as 


where  Xjj  are  the  components  of  the  stress  tensor,  are  the  components  of  the  strain 

tensor,  and  Cjjy  are  the  components  of  the  elasticity  tensor.  Note  that  Einstein’s  sum- 
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mation  convention  applies  to  equations  1.4  and  to  all  subsequent  equations  that  are 
written  in  tensor  notation,  unless  noted  otherwise. 


The  components  of  the  stress  and  strain  tensor  are  symmetric,  that  is 


^ij 


^kl  “  ^Ik’ 


(1.5) 


so  that 

‘^ijki  =  ^jiki  =  ‘^kiij  =  ^ijik- 

As  a  result,  there  are  at  most  36  independent  components  of  the  elasticity  tensor  in  the 
general  Hooke’s  law  for  anisotropic  materials.  The  existence  of  a  strain-energy  density 
function 


further  reduces  the  maximum  number  of  independent  components  to  21. 

Because  of  the  symmetries  in  the  components  Cjjki,  a  shortened  matrix  notation  has 
been  introduced  to  describe  the  stress-strain  relation  for  an  anisotropic  material  (Payton 
1983).  In  this  notation,  pairs  of  subscripts  involving  the  numbers  1, 2,  and  3  are  re¬ 
placed  by  1,  2,  3, 4,  5,  and  6  as  follows 

(11)  4-»l,  (22)  <^2,  (33)  4-»3,  (23)  =  (32)  4, 

(31)  =  (13)  4^  5  and,  (12)  =  (21)4-»6, 
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so  that,  for  example, 


^1232  “  ^64- 


(1.9) 


Using  the  shortened  notation,  we  write  equations  1.4  as 


= 


where  i,j  =  1,  2, ...,  6,  and 


11 

'*'22  ~  h 

■^33 

=  T 

'23 

II 

hi  ~  '^5 

^12 

=  T, 

and 


Ell  -  Ej 

^22  “  h  ^33  ~ 

2623  =  £4 

2631  =  £5  2£j2 

(1.10) 


(1.11) 


(1.12) 


The  components  Cjj  in  equations  1.10  are  called  the  elastic  moduli  of  the  material. 
When  the  elastic  moduli  are  independent  of  position,  the  medium  is  said  to  be  elastical¬ 
ly  homogeneous  and  the  components  Cjj  are  the  elastic  constants  of  the  material. 

When  the  material  is  elastically  symmetric  in  certain  directions,  the  number  of  in¬ 
dependent  elastic  moduli  in  equations  1 . 10  is  further  reduced.  For  instance,  when  there 
is  one  plane  of  elastic  symmetry  in  the  material,  there  are  13  independent  elastic  moduli 
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and  the  material  is  called  monoclinic.  When  there  is  elastic  symmetry  with  respect  to 
three  mutually  perpendicular  planes,  the  number  of  independent  elastic  moduli  reduces 
to  9  and  the  material  is  called  orthotropic.  If,  at  every  point  in  the  material,  there  is  one 
plane  in  which  the  elastic  moduli  are  equal  in  all  directions,  the  material  is  called  trans¬ 
versely  isotropic  and  the  number  of  elastic  moduli  reduces  to  5.  The  unit  vector  normal 
to  the  special  plane  of  isotropy  in  a  transversely  isotropic  material  is  called  the  preferred 
direction  of  the  material.  Comprehensive  discussions  of  material  symmetry  and  elastic 
moduli  can  be  found  in  Jones  (1975),  Christensen  (1979)  and  Sokolnikoff  (1987). 

1  ■2.2.2  Constraints  on  the  elastic  moduli 

The  elastic  moduli  Cy  are  positive  definite  when 

CijXiXj>0,  (i,j  =  1,2, ...,  6),  (1.13) 

for  arbitrary,  non-zero  vectors  Xj.  Physically  this  corresponds  to  the  requirement  that 
the  strain-energy  density  function  must  remain  positive  in  order  that  this  energy  be  min¬ 
imal  in  a  state  of  equilibrium.  The  necessary  and  sufficient  conditions  on  the  elements 
Cjj,  in  order  that  the  quadratic  form  (equation  1.13)  be  positive  definite,  are  that  all  of 
the  determinants. 


c 


11’ 


^11  ^12 
^21  ^22 


^11  ^12  ^13 
^21  ^22  ^^23 
^31  ^32  ^33 


(1.14) 
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remain  positive  (Kreyszig  1972).  The  specific  constraints  on  the  elastic  moduli  of  a 
transversely  isotropic  material  are  discussed  in  chapter  3. 

1. 2,2.3  Effective  elastic  constants  of  homogeneous  composite  materials 

When  we  speak  of  the  elastic  constants  of  a  homogeneous,  composite  material,  we 
really  mean  the  effective  elastic  constants  of  the  material.  Before  defining  the  effective 
elastic  constants  of  a  homogeneous,  composite  material,  it  is  first  necessary  to  explain 
what  is  meant  by  a  composite  material  versus  a  composite  body. 

A  composite  body  consists  of  two  or  more  materials  that  are  macroscopically  dis¬ 
tinct  and  are  joined  together  in  order  to  achieve  some  useful  property.  The  properties 
that  may  be  emphasized  include  wear  resistance,  thermal  insulation,  acoustic  insula¬ 
tion,  and  low  weight.  Examples  of  a  composite  body  include  bimetals,  clad  metals,  and 
laminated  glass. 

In  a  composite  material,  it  is  possible  to  define  a  representative  volume  element 
(Hashin  1983).  Representative  volume  elements  are  large  compared  to  the  dimensions 
of  the  components  of  the  composite  material.  A  necessary  characteristic  of  a  composite 
material  is  statistical  homogeneity.  In  a  statistically  homogeneous  composite  material 
all  global  characteristics  such  as  volume  fractions,  two-point  correlations,  etc.,  are  the 
same  in  any  representative  volume  element,  irrespective  of  its  position  within  the  ma¬ 
terial.  Examples  of  a  composite  material  include  fibrous  composites  and  particulate 
composites. 

The  effective  elastic  constants  of  a  composite  material  define  the  relations  between 
averages  of  field  variables  such  as  stress  and  strain,  when  the  spatial  variation  of  these 
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field  variables  is  statistically  homogeneous.  The  use  of  effective  elastic  constants  to  de¬ 
fine  the  relations  between  field  variables,  when  the  composite  material  is  subjected  to 
arbitrary  boundary  conditions  and,  therefore,  the  internal  field  variables  are  no  longer 
statistically  homogeneous,  is  discussed  by  Hashin  (1983)  and  is  summarized  in  the  fol¬ 
lowing  paragraphs. 

In  the  classical  continuum  mechanics  of  homogeneous  materials,  it  is  always  as¬ 
sumed  that  the  continua  retain  their  properties  regardless  of  specimen  size.  All  real  ma¬ 
terials,  however,  have  microstructure.  Metals,  for  example,  are  actually  polycrystalline 
aggregates  and  are  heterogeneous  materials  that  are  statistically  homogeneous.  There¬ 
fore,  the  differential  element  of  the  theory  of  elasticity  is  actually  a  representative  vol¬ 
ume  element.  The  effective  elastic  moduli  of  this  representative  volume  element  are  the 
elastic  moduli  of  the  theory  of  elasticity. 

The  classical  approach,  then,  in  the  analysis  of  homogeneous  composite  materials, 
is  to  assume  that  the  field  equations  of  elasticity  are  valid,  with  the  effective  properties 
of  the  composite  material  replacing  the  usual  elastic  properties  of  a  homogeneous  con¬ 
tinuum. 

1 .2.2.4  Wave  propagation  in  composite  materials 

The  dynamic  response  of  deformable  composite  materials  can  be  broadly  divided 
into  two  categories  (Christensen  1979).  If  the  wavelength  of  the  response  of  the  mate¬ 
rial  is  long  compared  to  the  scale  of  the  representative  volume  element,  the  material  re¬ 
sponse  is  governed  by  the  effective  properties  of  the  material.  In  this  case,  the  equations 
of  motion  of  the  composite  material  are  identical  to  those  of  homogeneous  materials. 
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If  the  wavelength  of  the  response  is  not  ideally  long  with  respect  to  the  representa¬ 
tive  volume  element,  very  complicated  dynamic  effects  occur.  These  effects  include 
wave  reflection  and  refraction  at  the  interfaces  of  the  components  of  the  composite  ma¬ 
terial.  These  phenomena  are  most  important  in  the  field  of  ultrasonics. 

In  this  study,  we  restrict  our  analysis  to  the  wavenumber  and  frequency  domains 
such  that  wavelengths  are  long  enough  so  that  we  can  assume  that  the  material  response 
is  governed  by  the  effective  material  properties. 

1.2.3  Previous  investigations  into  wave  propagation  in  cylinders 
1. 2.3.1  Shells  and  cylinders  in  a  vacuum 

Various  models  have  been  developed  to  investigate  the  free  wave  propagation  in 
hollow  circular  cylinders  in  a  vacuum.  Many  of  the  earlier  models  were  based  on  as¬ 
sumptions  that  restricted  their  validity  to  low  frequencies  of  shells  whose  thickness  is 
small  compared  to  the  inner  radius,  (Love  1944,  lunger  and  Rosato  1954,  Naghdi  and 
Berry  1954,  Lin  and  Morgan  1956a,  Naghdi  and  Cooper  1956,  Mirsky  and  Herrmann 
1957,  Cooper  and  Naghdi  1957,  Herrmann  and  Mirsky  1956).  Later,  Mirsky  and  Her¬ 
rmann  (1958)  investigated  wave  propagation  in  a  thick-walled  shell  but  restricted  their 
analysis  to  the  lowest  or  fundamental  mode. 

Other  approximate  models  have  been  developed  based  on  the  theory  of  elasticity. 
Bird  (1960)  and  Bird  et  al.  (1960)  used  a  perturbation  theory  that  lead  to  results  involv¬ 
ing  algebraic  and  trigonometric  functions.  McNiven  et  al.  (1966a)  expanded  the  dis¬ 
placements  in  a  series  of  orthogonal  polynomials.  McFadden  (1954)  and  Gazis  (1958) 
investigated  several  modes  of  the  plane  strain  vibrations  of  finite  or  infinite  cylinders 
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using  the  theory  of  elasticity,  but  excluded  displacements  in  the  axial  direction.  Finally, 
investigations  of  the  axisymmetric  and  nonaxisymmetric  modes  of  wave  propagation 
in  hollow  cylinders,  based  on  the  three-dimensional  equations  of  elasticity,  were  con¬ 
ducted  by  Armenakas  et  al.  (1969),  Gazis  (1959a,  b)  and  by  Kumar  and  Stephens 
(1972). 

Almost  all  of  the  above  investigations  were  restricted  to  the  purely  real  dispersion 
curves,  i.e.  involving  relationships  only  between  real  wavenumbers  and  real  frequen¬ 
cies.  However,  the  investigations  of  McNiven  et  al.  (1966b)  and  Kumar  and  Stephens 
(1972)  encompassed  purely  real,  purely  imaginary  and  complex  dispersion  curves. 

1. 2.3.2  Fluid-Filled  shells  and  cylinders 

Wave  propagation  in  fluid-filled  tubes  has  been  investigated  since  the  time  of 
Young  (1808),  who  produced  a  comprehensive  analysis  of  the  influence  of  wall  com¬ 
pliance  on  blood  flow  in  human  arteries.  Since  that  time,  the  subject  of  blood  flow  has 
been  studied  extensively,  (see  Skalak  (1966)  for  a  survey  of  the  literature  on  blood 
flow).  Examples  of  other  studies  includes  that  of  Fay  et  al.  (1947)  and  Jacobi  (1949) 
who  used  various  approximations,  such  as  a  restriction  to  low  frequencies  or  ignoring 
the  Poisson’s  coupling  of  the  shell,  in  their  investigations  of  axisymmetric  wave  prop¬ 
agation  in  tubes.  Thomson  (1953)  introduced  the  effects  of  Poisson’s  ratio  and  evaluat¬ 
ed  the  phase  velocities  of  the  first  three  axisymmetric  “fluid”  waves.  Later,  Lin  and 
Morgan  (1956b)  studied  the  axisymmetric  wave  propagation  in  a  fluid-filled  shell  with¬ 
out  the  approximations  introduced  by  previous  authors. 

Kumar  (1971, 1972)  investigated  the  axisymmetric  and  nonaxisymmetric  wave 
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propagation  in  a  fluid-filled  cylinder  of  arbitrary  thickness,  within  the  framework  of  the 
theory  of  elasticity,  and  analyzed  the  cut-off  frequencies  and  the  real,  imaginary,  and 
complex  dispersion  curves  as  functions  of  wall  thickness.  More  recently.  Fuller  and 
Fahy  (1982)  investigated  the  real,  imaginary,  and  complex  dispersion  behavior  and  en¬ 
ergy  distributions  of  axisymmetric  and  nonaxisymmetric  waves  in  thin-walled  cylindri¬ 
cal  elastic  shells,  filled  with  a  fluid,  as  functions  of  frequency  and  material  parameters. 
1. 2.3.3  Shells  and  cylinders  immersed  in  a  fluid 

Many  of  the  early  investigators  of  the  vibrations  of  thin  and  thick  shells  surrounded 
by  a  fluid  were  interested  in  the  application  of  the  results  to  acoustic  waveguides  and 
underwater  sound  sources.  lunger’s  (1952a,  b)  early  analysis  of  thin  shells  considered 
only  axisymmetric,  radial  motion  of  the  shell.  lunger  (1953)  later  considered  the  three- 
dimensional  thin  shell  problem.  In  this  study,  lunger  did  not  address  the  fluid-solid  in¬ 
teraction  problem  but,  instead,  assumed  that  the  shell  had  known  displacements.  Bleich 
and  Baron  (1954)  extended  the  work  of  lunger  (1953)  by  utilizing  the  modes  of  vibra¬ 
tion  of  the  thin  shell  in  a  vacuum  as  generalized  coordinates  and  emphasized  the  struc¬ 
tural  response  of  the  shell  rather  than  the  acoustic  response  of  the  fluid.  Greenspon 
(1960)  studied  the  case  of  the  axisymmetric  vibration  of  a  thick  cylindrical  shell  in  wa¬ 
ter  for  the  purpose  of  gaining  insight  into  the  shell  behavior  and  specialized  the  results 
for  applications  to  underwater  sound  sources  and  transducers.  Greenspon  (1961)  ex¬ 
tended  this  investigation  into  the  more  general  type  of  deformation  corresponding  to 
higher  order  circumferential  modes  of  vibration  of  thin  and  thick  shells  in  water.  The 
text  by  lunger  and  Feit  (1986)  provides  a  comprehensive  study  on  the  vibration  of 
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plates  and  shells  in  a  vacuum  and  in  water  and  of  the  sound  radiation  into  the  fluid. 

1 .2. 3 .4  Shells  and  cylinders  that  are  fluid-filled  and  immersed  in  a  fluid 

Due  to  the  complexity  of  the  geometry,  there  have  been  few  investigations  into  the 
wave  propagation  in  circular  cylinders,  of  arbitrary  thickness,  with  fluid  loading  on  the 
inner  and  the  outer  surfaces.  Chandra  and  Kumar  (1976)  investigated  the  dispersion  of 
axisymmetric  waves  in  fluid-filled  cylinders  that  had  external  fluid  loading.  They  found 
that  the  effect  of  the  surrounding  fluid  on  the  fluid-filled  cylinder  increased  with  a  de¬ 
crease  in  the  wall  thickness  of  the  cylinder  and  extended  their  results  into  the  complex 
wavenumber  domain.  Sinha  et  al.  (1992)  also  investigated  the  case  of  axisymmetric 
wave  propagation  in  a  fluid-loaded  cylinder  of  fixed  wall  thickness.  In  their  study,  the 
boundary  conditions  increased  in  complexity  from  a  hollow  cylinder  in  a  vacuum  to  a 
fluid-filled  cylinder,  a  hollow  cylinder  with  fluid  loading  on  the  outer  surface,  and  a 
cylinder  with  fluid  loading  on  the  inner  and  the  outer  surfaces.  A  comparison  of  the  an¬ 
alytical  results  were  made  to  experimental  data  by  Fiona  et  al.  (1992). 

1. 2.3.5  Anisotropic  cylinders 

All  of  the  investigations,  discussed  so  far,  have  dealt  with  cylinders  composed  of 
isotropic  materials.  A  limited  amount  of  investigations  exists  in  the  area  of  wave  prop¬ 
agation  in  cylinders  composed  of  anisotropic  materials,  particularly  based  on  the  three- 
dimensional  theory  of  elasticity.  Chree  (1890)  may  have  been  the  first  investigator  of 
axisymmetric  wave  propagation  in  an  anisotropic  bar.  In  this  analysis,  the  equations  of 
motion  of  a  transversely  isotropic  bar,  of  any  cross  section,  were  developed  and  then 
specialized  for  the  case  of  a  circular  rod.  Approximate  solutions  were  obtained  by  ex- 
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ponding  the  displacements  in  powers  of  the  radius  r.  Later,  Morse  (1954)  developed  the 
exact  solution  of  the  problem  studied  by  Chree  (1890)  and  showed  that  the  solution  re¬ 
duced  to  the  Pochhammer  solution  for  the  isotropic  case.  No  numerical  results  were 
presented  by  either  Chree  (1890)  or  Morse  (1954).  Mirsky  (1965a,  b)  investigated  axi- 
symmetric  and  nonaxisymmetric  wave  propagation  in  transversely  isotropic  cylinders 
and  used  displacement  potentials  to  solve  the  equations  of  motion.  Fraser  (1980)  used 
the  method  of  eigenfunction  expansion  to  separate  the  equations  of  motion  of  an  aniso¬ 
tropic  cylinder  and  specialized  the  results  to  a  transversely  isotropic  material.  However, 
he  did  not  analyze  the  dispersion  of  the  modes  of  vibration.  Axisymmetric  wave  prop¬ 
agation  in  orthotropic  cylinders  was  studied  by  Mirsky  (1964)  who  obtained  closed- 
form  solutions  of  the  equations  of  motion  in  terms  of  infinite  series  using  the  method 
of  Frobenius.  Armenakas  and  Reitz  (1973)  expanded  the  displacements  of  the  cylinder 
into  a  power  series  of  the  radial  coordinate  to  study  wave  propagation  with  an  arbitrary 
number  of  circumferential  modes  in  an  orthotropic  cylinder.  However,  none  of  the 
above  investigations  included  the  effects  of  fluid  loading  on  either  the  internal  or  exter¬ 
nal  surfaces  of  the  cylinder. 

To  the  author’s  knowledge,  wave  propagation  in  fluid-loaded,  transversely  isotro¬ 
pic  cylinders  has  not  been  investigated,  to  date,  within  the  framework  of  the  theory  of 
elasticity, 

1.3  OVERVIEW  OF  THE  DISSERTATION 

In  chapter  2,  the  displacement  equations  of  motion  of  the  cylinder  are  formulated 
using  the  constitutive  equations  of  a  transversely  isotropic  material  with  a  preferred  ma- 
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terial  direction  collinear  with  the  longitudinal  axis  of  the  cylinder.  The  displacement 
equations  of  motion  of  the  internal  and  external  fluids  are  formulated  using  the  consti¬ 
tutive  equations  of  an  irrotational,  inviscid  fluid.  Displacement  potentials  are  used  to 
solve  the  equations  of  motion  of  the  cylinder  and  of  the  fluids.  The  frequency  equation 
of  the  coupled  system,  consisting  of  the  cylinder  and  the  inner  and  outer  fluids,  is  de¬ 
veloped  under  the  assumption  of  perfect-slip  boundary  conditions  at  the  fluid-solid  in¬ 
terfaces.  This  frequency  equation  is  general  in  axial  wavenumber  k,  circumferential 
wavenumber  n,  wall  thickness  h,  and  radian  frequency  to. 

In  chapter  3,  the  frequency  equation  is  specialized  for  the  case  of  material  isotropy 
and  for  the  n=l,  nonaxisymmetric  modes  of  wave  propagation.  Cut-off  frequencies  are 
computed  for  the  first  four  plane  strain  modes  and  the  first  two  axial  shear  modes,  as 
functions  of  two  different  Poisson’s  ratios  and  the  ratio  of  the  cylinder  inner  to  outer 
radius.  In  particular,  Poissons ’s  ratios  for  steel  and  for  rubber  are  used  in  these  compu¬ 
tations.  The  frequency  spectrum  of  wave  propagation  in  a  hollow  cylinder,  a  hypothet¬ 
ical  fluid  column,  a  fluid-filled  cylinder,  and  a  cylinder  that  is  fluid-filled  and  immersed 
in  fluid,  is  computed  for  each  of  the  two  materials. 

In  chapter  4,  the  cut-off  frequencies  and  frequency  spectrum  are  computed  for  the 
n=l,  nonaxisymmetric  modes  of  wave  propagation  in  a  homogeneous,  transversely  iso¬ 
tropic  cylinder  that  is  fluid-filled  and  immersed  in  fluid.  The  particular  geometry  that 
is  investigated  is  a  fiber-reinforced  cylinder  with  the  fibers  aligned  parallel  to  the  lon¬ 
gitudinal  axis  of  the  cylinder.  The  effective  elastic  constants  are  computed  using  the  in¬ 
dividual  elastic  constants  of  the  matrix  and  reinforcing  fibers  and  the  volume  fractions 


15 


of  each  of  the  two  components.  The  numerical  results  are  compared  to  those  obtained 
for  the  two  isotropic  materials. 

A  summary  and  conclusions  of  the  results  of  chapters  2  to  4  are  provided  in  chapter 
5  and  a  discussion  of  recommendations  for  continued  investigations  are  provided  in 
chapter  6. 

The  Mathematica  programming  language  (Wolfram  1991)  was  used  extensively 
throughout  this  investigation.  The  symbolic  computation  capabilities  of  Mathematica 
were  first  used  to  assist  in  the  manipulation  of  algebraic  formulas  in  the  analytical  for¬ 
mulation  of  the  problem.  Next,  Mathematica’s  numerical  capabilities,  including  com¬ 
putations  of  complex  Bessel  functions,  were  used  for  the  root  searches  of  the  complex 
fr^uency  equation. 
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2.  ANALYTICAL  FORMULATION 

We  begin  this  chapter  with  the  physical  description  and  assumptions  of  the  behavior 
of  the  coupled  system  consisting  of  a  hollow  cylinder  filled  with  and  immersed  in  an 
inviscid  fluid.  We  then  develop  the  displacement  equations  of  motion  of  the  cylinder 
using  the  constitutive  equations  of  a  linearly  elastic,  homogeneous,  transversely  isotro¬ 
pic  material,  and  the  assumed  displacement  forms.  Solutions  of  the  equations  of  motion 
are  used  to  derive  the  displacement  components  and  surface  stresses  on  the  cylinder. 
Next,  we  develop  the  displacement  equations  of  motion  of  the  internal  and  external  flu¬ 
ids,  using  the  constitutive  equations  of  an  irrotational,  inviscid  fluid,  and  the  assumed 
displacement  forms.  Solutions  of  the  equations  of  motion  of  the  internal  and  the  exter¬ 
nal  fluids  provides  their  respective  displacement  components  and  pressures. 

The  frequency  equation  of  the  system  is  developed  by  coupling  the  cylinder  to  the 
fluids  through  the  boundary  conditions  at  the  cylinder’s  inner  and  outer  surfaces.  The 
classical  approach  is  used  to  deal  with  the  inviscid  fluid-solid  interface.  That  is,  perfect- 
slip  boundary  conditions  are  assumed  allowing  for  discontinuity  in  the  planar  displace¬ 
ment  components.  Application  of  these  boundary  conditions  results  in  a  system  of 
eight,  linear,  homogeneous,  algebraic  equations  in  the  unknown  wave  amplitudes.  The 
solutions  of  these  equations  are  discussed  in  chapters  3  and  4. 

2.1  SYSTEM  GEOMETRY 

The  system  consists  of  an  infinitely  long,  hollow,  linearly  elastic,  transversely  iso¬ 
tropic  cylinder  with  inner  radius  a,  and  outer  radius  b.  The  cylinder  is  filled  with  an  ir- 
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rotational,  inviscid  fluid  with  a  density  and  an  acoustic  phase  velocity  c^ .  The 

f 

cylinder  is  also  immersed  in  a  second,  irrotational,  inviscid  fluid  with  a  density  P2  and 
an  acoustic  phase  velocity  C2 .  The  system  is  assumed  to  be  linear  so  that  the  linearized, 
three-dimensional,  stress  equations  of  motion  can  be  utihzed  for  the  cylinder  and  the 
fluids.  The  system  displacements  and  stresses  are  defined  by  the  cylindrical  coordi¬ 
nates,  r,  8,  z.  The  cylindrical  coordinate  directions  are  designated  as  the  1,  2,  and  3  di¬ 
rections,  respectively. 

2.2  EQUATIONS  OF  MOTION  OF  THE  CYLINDER 
2.2.1  Three-dimensional  dynamic  equations  of  elasticity 

The  three-dimensional  stress  equations  of  motion  and  strain-displacement  relation¬ 
ships  of  a  linearly  elastic,  solid  medium  are 


=  pu’i- 


(2.1) 


and 


=  -  (u-  •  +  u- .) 
u  2  ^ 


(2.2) 


respectively,  where  Ty  are  the  components  of  the  stress  tensor,  are  the  components 
of  the  strain  tensor,  Uj  are  the  displacement  components,  and  p  is  the  density.  It  should 
be  noted  that  Einstein’s  summation  convention  applies  to  equations  2.1  and  2.2  and  to 
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all  other  equations  that  are  written  in  tensor  notation,  and  that  Xjj  j  denotes  3Xjj/3xj , 


etc. 


Equations  2.1  and  2.2  are  written  in  cylindrical  coordinates  as 


Vr  +  7Vo  +  Vz  +  ^''r,-W  =  P«, 

'^ez,z  j ^re  “^re.r  “ 

Vz  +  ;Ve  +  H^  + Vr  =  P“z’ 


(2.3) 


and 


e 


r 


U 


r,r 


1 

^0  “  j.  ■'■7^0.6 


e 


Z 


U 


z,z 


2e,e  =  2ee,  =  Ue^-^  +  ^u^e 
2e0z  =  2e,e  =  ^u^.e  +  ug,. 


2erz  = 


u  +  u 
r,z  z,r 


(2.4) 


respectively. 

2.2.2  Constitutive  equations 

The  constitutive  equations  of  a  general  anisotropic  material  are 


(2.5) 
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where  Cjjy  are  the  components  of  the  elasticity  tensor.  As  discussed  in  chapter  1,  when 
the  medium  is  homogeneous  the  components  of  the  elasticity  tensor  are  constant  and, 
when  the  medium  is  homogeneous  and  transversely  isotropic  only  5  of  the  elastic  con¬ 
stants  are  independent. 

Homogeneous,  transversely  isotropic  materials  have  one  plane  in  which  the  elastic 
constants  are  equal  in  all  directions.  The  unit  vector,  normal  to  the  special  plane  of  isot¬ 
ropy  in  a  transversely  isotropic  material,  is  called  the  preferred  direction  of  the  material. 
When  the  preferred  direction  of  a  transversely  isotropic  material  is  given  by  a  unit  vec¬ 
tor  {m},  the  constitutive  equations,  2.5,  can  be  written  as  (Sutcu  1992) 


"ij 


=  Sij(Ciekk  +  C2mpmqepq)  +  +C^e.. 


(2.6) 


where  the  constants  Cj,  C2,  C3,  C4,  and  C5,  are  related  to  the  more  familiar  engineering 
constants  as 

Cl  =  By  —  Gj 

C2  =  (2Vl-  l)By  +  Gy 

C3  =  El+  (2vL-l)^By-4GL  +  Gy  (2.7) 

C4  =  2Gy 

C5  =  2GL-2Gy. 

The  definitions  of  the  engineering  constants  on  the  right-hand  side  of  equations  2.7  are 
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=  the  transverse  bulk  modulus, 

Gl  =  the  longitudinal  shear  modulus, 

Gx  =  the  transverse  shear  modulus, 

El  =  the  longitudinal  Young’s  modulus. 

Ex  =  the  transverse  Young’s  modulus, 

Vl  =  the  longitudinal  Poisson’s  ratio, 

Vx  =  the  transverse  Poisson’s  ratio, 

where  “transverse”  refers  to  properties  in  the  plane  of  isotropy  and  “longitudinal  refers 
to  properties  along  the  material  preferred  direction. 

When  the  preferred  direction  vector  {m}  is  collinear  with  the  z,  or  the  3,  direction 
of  the  cylinder,  i.e.  mi=0,  m2=0,  and  m3=l,  equations  2.6  become 


where 


'^rr  “  ^12^00 


'^00  “  ^12^rr'’’‘^11^00  ^13^zz 


(2.8) 


'^zz  ”  ^13^00  "^^SS^zz 


''0Z  “ 


2^44^0z 


'^zr  ~  ^^44^zr  '^'10 


“  2^66^10  ’ 


'^11  “ 


'^22  ~  “^00 


'^33  “  "^zz 


■^31  “  '^zi 


'^23  ~  ^0z 


'^12  ~  '^'10  ’ 


(2.9) 


^11 

=  Err 

^22 

~  ^ee 

^33  -  ^zz 

^23 

=  Eez 

£31 

=  ^zr 

^12  “  ^r0  ’ 

and  where 


'11 


'12 


'13 


=  Cj  +  =  B-p  +  Gj 

=  Cj  =  Bj  — Gq^ 

=  Cj  +  C2  =  2BjVp^ 


c 

c 


23  =  Cj  +  2C2  +  C3  +  +  2C5 


(C4  +  C5) 


44 


=  G. 


=  El  +  4B-pVl 


'66 


-  —  Gj  -  (Cjj  Cj2)/2. 


(2.10) 


(2.11) 


A  Mathematica  program  (Wolfram  1991)  was  developed  to  expand  equations  2.6 
for  any  unit  vector  {m}  and  to  provide  the  resulting  elastic  constants  in  terms  of  the  en¬ 
gineering  constants.  This  program  is  provided  in  appendix  D. 

Substitution  of  the  constitutive  equations,  2.8,  and  the  strain-displacement  relation¬ 
ships,  2.4,  into  the  stress  equations  of  motion,  2.3,  results  in  the  following  displacement 
equations  of  motion. 
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'11 


^1 

'^r.rr  2 

V  ^  r  y 


1 


^66  .1 


T“r.ee  ^44*^1,22  j  (^66  ^12)  “e,re 


2  ^^66  '*’  ^11^  ^0,0  '*’  ^^44  "^^13^  '^z.rz  P^r 


;  (<=66  +  C12)  “,.,e  +  1  (<=66  +  <=ll)  “,.6  +  ‘=66 
^  r 

+  ^UQ00  +  C44Ue,^  +  J(c44  +  Ci3)u^e,  =  piig 


Ue.rr  +  ;Ue^r-^ 

r  y 


(C44  +  C13)  u^  ^  (C44  +  C13)  ^  (C44  +  c,3)  Ue  0^ 

+  C44(u,.„+^U,,,  +  4«z.ee)  +  C33'^2.zz  =  PV 


(2.12) 


2.2.3  Displacement  potentials 

To  uncouple  equations  2.12,  we  follow  Mirsky  (1965)  who  assumed  the  following 
displacement  forms, 


u^  (r,  e,  z,  t)  =  [<[)  (r,  0) ,,  +  (r,  6)  ,0  jcos  (cot  +  kz) 

U0  (r,  0,  z,  t)  =  Q<|>  (r,  0)  ,e  -  V  (r,  0)  .^jcos  (cot  +  kz) 

u^  (r,  0,  z,  t)  =  Ti<t)  (r,  0)  sin  (cot  +  kz) , 

where  (|)(r,0)  and  \|i(r,0)  are  displacement  potentials  and  T)  is  an  arbitrary  constant. 
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For  mathematical  convenience  we  will  use  an  exponential  dependence  on  z  and  on 
t  so  that  the  assumed  displacement  forms  become 


Uj(r,0,  z,t)  =  (^(t)(r,0),^  +  |\|;(r,  e),Qje 

ue(r,0,z,t)  =  Q<t)(r,0),e-v(r,0),^y*^“‘'’'"^  (2-14) 

u^  (r,  0,  z,  t)  =  r\^  (r,  0)  . 


With  this  notation,  complex  phase  velocities  can  be  handled  in  a  straightforward  man¬ 
ner. 

When  equations  2.14  are  substituted  into  equations  2.12,  and  the  results  are  simpli¬ 
fied,  the  equations  of  motion  of  the  displacement  potentials  are  found  to  be 


[c„V2<t)(r,0)  +  {pco^-k^C44  +  iTik(Cj3-i-C44)}(})(r,  0)j,^ 
+  “[c66^V(r.0)  +(pa)^-k^C44)\|/(r,0)],e  =  0 


^j^CjiV2(t)(r,0)  -I-  {po)^-k^C44-HiTik(Ci3  +  C44)}(|)(r,0)j,e 

-[c66^V(r.0)  +(pa)^-k^C44]\|/(r,0)],^  =  0 

[ik  (Ci3  +  C44)  +  TIC44]  V2(t)  (r,  0)  -I-  po)^  -  k%3]Ti(t)  (r,  0)  =  0, 


(2.15) 
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where 


v"*  =  (  ).„  +  J(  ).,+  -2(  ).66  (2.16) 

^  r 

is  the  two-dimensional  Laplacian  operator,  and  where  the  exponential  dependence  on  t 
and  z  has  been  suppressed. 

The  first  two  equations  of  the  set  2.15  are  satisfied  if,  in  particular, 

CjiV2(t)(r,  0)  +  [pto^-k\44  +  irik(Cj3  4-C44)](>(r,  0)  =0,  (2.17) 


and  if 


CggV^YCr,  0) +(^pa)^-k^C44j\i/(r,  0)  =0.  (2.18) 

The  last  equation  of  the  set  2.15  is  consistent  with  equation  2.17  when  Ti  is  the  solution 
of 


Tljpco^-kSg] 

ik(Ci3-HC44)  +110^ 


2  2 

pco  -k  C44-l-iTlk(Ci3 -I-C44) 


'11 


(2.19) 


After  solving  for  H  from 


25 


(2.20) 


pco^  -  k^C44  +  iTik  (Ci3  +  C44) 


'11 


and  substituting  the  solution  into 


Ti[po)^-k%3j  ^  2 

ik(Cj^3+C^)  + 

we  find  that  ^  satisfies  the  equation 


(2.21) 


(CiiC^^)^  [(Cij+C^)pCO  +  (^Cj3  +  2CJ3C44  ^11^33 


+  ipO)^-C33k^]  pC0^-C44k^J  =  0. 


(2.22) 


The  solutions  of  are 


c 


2 

2 


-B-Vb^-4AC 

2A 

-B  + Vb^-4AC 
2A 


(2.23) 


where 
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(2.24) 


From  equations  2.17,  the  third  equation  of  the  set  2.15  and  equation  2.18,  the  equa¬ 
tions  of  motion  of  the  displacement  potentials  (|)i(r,e)  (corresponding  to  from  equa¬ 
tion  2.20),  <t>2(r,6)  (coiresponding  to  from  equation  2.21),  and  \j/(r,0),  can  be  shown 
to  be 


V2<t)^(r,e)  +c?<l>i(r,e)  =  0 

V2(t)2(r,e)+C2<l>2(r^e)  =0 

V2v(r,0)  +qV(r,e)  =  0. 


where 


2 

q 


'66 


(2.25) 


(2.26) 
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2.2.4  Solutions  of  the  equations  of  motion  of  the  cylinder 

We  first  consider  the  solution  of  <!)j  (r,  0) .  In  view  of  the  finite  cross-sectional  di¬ 
mensions  of  the  cylinder,  we  assume  a  separation  of  variables  solution  in  r  and  0 ,  i.e. 

4)1  (r,0)  =  Oi(r)0i(0).  (2.27) 

Substitution  of  equation  2.27  into  the  first  equation  of  set  2.25,  results  in 


(r) , 


rr 


0i(0),ee  +  n^0i(0)  =  0 

2 

^  r 


(2.28) 


The  solutions  of  the  first  equation  in  set  2.28  are  sines  and  cosines  of  argument  n0 . 
Because  the  solutions  should  be  continuous  functions  of  0  with  continuous  derivatives, 
n  can  only  be  zero  or  an  integer.  The  second  equation  of  set  2.28  is  Bessel’s  equation 
of  order  n  when  is  real  or  complex,  and  is  the  modified  Bessel’s  equation  of  order  n 
when  ^1  is  imaginary.  The  solution  of  Oj  (r)  is 

®i(r)  =  A,W„(4,r) +B,Z„(5,r),  (2.29) 
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where  =  Ci  and  and  are  the  Bessel  functions  Jn  and  ¥„,  respectively,  when 
is  real  or  complex,  and  =  |Ci|  ^  ^  modified  Bessel  functions 

Iji  and  Kji,  respectively,  when  is  imaginary. 

Assuming  a  separation  of  variables  solution  for  (])2(r,0)  and  V|/(r,6)  in  r  and  0  results 

in 


<I>2  (r)  =  AjW,  (^jr)  +  B2Z„  (IjI)  ,  (2.30) 

and 

T  (r)  =  A3W„  (y)  +  B3Z„  i^) .  (2.31) 

In  equation  2.30,  =  C2  and  Z^  are  the  Bessel  functions  Jn  and  Yn,  respec¬ 
tively,  when  ^2  is  real  or  complex,  and  ^2  =  IC2I  ^  modified 

Bessel  functions  and  Kn,  respectively,  when  ^2  is  imaginary.  In  equation  2.31, 

^3  =  q  and  Wn  and  Zn  are  the  Bessel  functions  Jn  and  Yn,  respectively,  when  q  is  real 
or  complex,  and  ^3  =  |q|  and  Wn  and  Zn  are  the  Bessel  functions  In  and  Kn,  respec¬ 
tively,  when  q  is  imaginary. 

The  sine  or  cosine  dependence  on  0  in  (t)i(r,0),  (l)2(r,0),  and  \ir(r,0)  must  be  chosen 
so  that  they  satisfy  equation  2.14.  It  is  evident  that  a  cosine  dependence  on  0  in  (t)i(r,0) 
and  in  (t)2(r,0)  is  consistent  with  a  sine  dependence  on  0  in  \|r(r,0),  and  vice  versa.  In  this 
investigation,  we  use  the  following  set  of  solutions  of  the  displacement  potentials. 
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(t)^(r,e)  =  (AiW„(^jr)+BiZJ^ir))cos(ne), 


(2.32) 


<1)2  (r,  0)  =  (A2W„  i^^t)  +  a^T) )  cos  (nB) ,  (2.33) 

and 

V  (r,  0)  =  (A3W„  (^jt)  +  (^3r) )  sin  (n0) .  (2.34) 

2.2.5  Cylinder  displacements  and  surface  stresses 

When  equations  2.32,  2.33  and  2.34  are  inserted  into  equations  2.14,  the  cylinder 
displacements  become 


u,(r,0,z,t)  =  {AiW„(^ir)  +BiZ„(^ir) 

**  *r 

+  ;[A3W„(l3r)  +B3Z„(43r)]  }  cos  (nG) 


(2.35) 
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(2.36) 


Ue(r,e,z,t)  =  {-2[A,W„(4ir) +B,Z„(|,r)] 
-;IA2W„(5jr)+B2Z,(42r)l 
-  [  A3W„  (^jr)  ^  +  BjZ^  (^jO  }  sin  (nS) 


and 


u^(r,0,z,t)  =  {Tij[AiW„(^ir) +BiZ„(^jr)] 
+  TI2 [ A2W„ (^2^)  +B2Zj,(^2^)]  }  cos(n0)e"‘^“‘ 


(2.37) 


where,  from  equations  2.20  and  2.21,  we  have 


Til  = 


TI2  = 


IkCi  (^13  ^44) 

€33^-044^1) 

f  2  2  2 

il-C44k  +p0)  -C11C2 


k 


(2.38) 
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Next,  when  equations  2.35, 2.36,  and  131  are  inserted  into  equations  2.4  and  2.8, 


and  simplified,  the  surface  stresses,  x„,  and  Xre  become 


=  {2c«;(A3W.(^3r)^  +  B3Z„(53r>.,-;[A3W„(V) 

+  B3Z„(^3r)  ]  )-2?^[A,W„a,r)  ,  +  B,Z,(^,r)^ 

2 

+  ,  +  B3Z„(4jr)  +2c„!^  [A.W„(^,r)  +B,Z„(|,r) 

’  ’  r 

+  BiZn(^ir)  ]  +  (iCi3'n2k-CjjC2) 


cos  (n0)e 


A2W„(^2r)  +B2Z„(^2^)  ]  } 

-i((ot-kz) 


=  c„{(Tl,+ik)  [A,W„(V),,  +  B,Z„(^,r)^] 
+  (ll3  +  ik)  [AjW„(42r)^  +  B3ZJ^3r)  J 
+  ik![A3W„(^3r)  +B3Z„(^3r)]  }  cos  (nB) e'*'”'" 


and 


(2.39) 


(2.40) 
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\e  =  <=66{2-2[A,W„(4,r)  +B,Z„(5,r)  ] 

r 

-2"  [A,W„  (^,r)  ,  +  BiZ„  (lir)  +  A^Wi^  (^ji)  +  B^Z„  (i^r)  ^  ] 

+  q'[A3W„  i^)  +  B3Z„  (^3!)  1  +  ?  |A3W„  a,t)  ^  +  BjZ^  a,r)  ,  1 

2 

-222|^A3W„(^3r)  +B3Z^(^3r)  ]  }sin(ne)e 


(2.41) 


-i(cot-kz) 


2.3  EQUATIONS  OF  MOTION  OF  THE  FLUID 
2.3.1  Stress  equations  of  motion  and  acoustic  pressure 

The  three-dimensional  stress  equations  of  motion  of  and  excess  acoustic  pressure 
in  a  linear,  irrotational,  inviscid  fluid  are 


(2.42) 


and 


P  = 


(2.43) 


respectively. 

2.3.2  Constitutive  equations 

The  constitutive  equations  of  an  inviscid  fluid  are 
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(2.44) 


=  -pS,, 

We  note  that  the  stress  tensor  for  an  inviscid  fluid  is  isotropic.  That  is,  the  components 
of  the  stress  tensor,  in  any  rectangular  Cartesian  system,  are  unaltered  by  orthogonal 
transformations  of  the  coordinates. 

Substitution  of  equations  2.43  and  2.44  into  2.42  provides  the  displacement  equa¬ 
tion  of  motion 


B'u'ji  =  p'u[,  (2.45) 

where  is  the  adiabatic  bulk  modulus  of  the  fluid. 

2.3.3  Displacement  potentials 

In  investigations  of  wave  propagation  in  an  inviscid  fluid,  it  is  generally  assumed 
(Kumar  1972,  Chandra  and  Kumar  1976,  Sinha  et  al.  1992)  the  displacements  to  be  of 
the  form. 


u[  (r,  e,  z,  t)  =  d/  (r,  0,  z,  t)  ^  (2.46) 

where  it  should  be  noted  that,  because  the  fluid  is  inotational,  (t)^(r,0,z,t)  must  be  a  scalar 
function. 

Upon  substitution  of  2.46  into  2.45,  it  can  be  shown  that  (t)^(r,0,z,t)  must  satisfy  the 
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wave  equation,  that  is 


(r,  0,  z,  t)  0  =  ® ' 

c 


(2.47) 


where 


2 

c  = 


f 

P 


(2.48) 


is  the  acoustic  phase  velocity  in  the  fluid,  and  where,  in  cylindrical  coordinates. 


(2.49) 


2.3.4  Solution  of  the  equations  of  motion  of  the  fluids 

The  radial  displacement  of  the  fluids  and  of  the  cylinder  must  be  continuous  at  the 
fluid-solid  interfaces.  Upon  inspection  of  equations  2.35  and  2.46,  it  is  evident  that 
<}>^(r,0,z,t)  must  be  of  the  form 
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(2.50) 


(t)^(r,  0,  z,  t)  =  0^(r)  cos  (n0)  e 


When  equation  2.50  is  substituted  into  2.47,  we  have 


‘'’'('•).rr  +  7*'W.r  + 


(0  ,2n^ 

-TT-k - r 

2  2 
VC  T  J 


d>^(r) 


cos (n0) e 


-i  (cot-kz) 


=  0.  (2.51) 


Equation  2.51  must  be  valid  for  all  0,  t,  and  z,  so  that 


^  2 

-1 

<0  , 2  n-^ 

— pr-k  — r 

2  2 

0^(r) 

V  c 

T  J 

=  0. 


(2.52) 


2  2  2 

Equation  2.52  is  Bessel’s  equation  of  order  n  when  to  /c  -  k  is  real  or  complex,  and 

2  2  2 

is  the  modified  Bessel’s  equation  of  order  n  when  to  /c  -  k  is  imaginary. 

2.3.4. 1  Inner  fluid 

The  solution  of  equation  2.52  for  a  fluid  column  must  remain  finite  as  r  0 .  The 
Bessel  function  of  the  second  kind  (^^r) ,  and  the  modified  Bessel  function  of  the 
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second  kind  Kjj  (^4^)  ,  are  unbounded  when  r  =  0 .  Therefore,  the  solution  of  equation 
2.52  for  the  inner  fluid  is 


=  DjWj44r),  (2-53) 

where  ^4  =  ®  and  Wn  is  the  Bessel  function  when  a  is  real  or  complex,  and 
^4  =  |ot|  and  Wji  is  the  modified  Bessel  function  when  a  is  imaginary,  and  where 


2 

a  = 


(2.54) 


2.3.4.2  Outer  fluid 

The  solution  of  equation  2.52  for  outgoing  waves  in  an  infinite  fluid  must  remain 
finite  as  r  00 .  The  modified  Bessel  function  (^jr)  is  unbounded  when  r  00 . 
Therefore,  the  solution  of  equation  2.52  for  the  outer  fluid  is 

02  (r)  =  D2W„(^5r),  (2.55) 

where  ^5  =  P  Wn  is  the  Hankel  function  of  the  first  kind  ,  when  p  is  real  or 
complex,  and  ^5  =  |p|  and  is  the  modified  Bessel  function  K^,  when  P  is  imagi¬ 
nary,  and  where 
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(2.56) 


2.3.5  Displacements  and  pressures  in  the  fluids 

The  relevant  displacement  and  pressure  components  for  the  inner  and  outer  fluids 
are  as  follows. 

2.3.5. 1  Inner  fluid 

When  equations  2.50  and  2.53  are  substituted  into  equation  2.46,  the  radial  compo¬ 
nent  of  displacement  for  the  inner  fluid  becomes 


u^^(r,e,z,t)  =  D^W^C^^r)  cos  (nO) (2.57) 

and,  when  equations  2.44  and  2.57  are  substituted  into  equation  2.42,  the  pressure  in 
the  inner  fluid  becomes 

Pi  =  p;o)^DiW„(^4r)cos(ne)e"‘^“‘“‘^"\  (2.58) 

2.3.5.2  Outer  fluid 

When  equations  2.50  and  2.55  are  substituted  into  equation  2.46,  the  radial  compo¬ 
nent  of  displacement  for  the  outer  fluid  becomes 
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(2.59) 


u[  (r,  0,  z,  t)  =  D2W„  i^^T)  ^cos  (n0)  e 

and,  when  equations  2.44  and  2.59  are  substituted  into  equation  2.42,  the  pressure  in 
the  outer  fluid  becomes 

P2  =  cos  (n0)e  .  (2.60) 

2.4  BOUNDARY  CONDITIONS 

To  complete  the  analytical  formulation  of  the  frequency  equation  of  the  coupled 
system  requires  consideration  of  the  boundary  conditions  at  the  inner  and  outer  surfaces 
of  the  cylinder.  In  a  solid-solid  interface  problem,  the  continuity  conditions  require  that 
the  three  displacement  components  and  the  three  surface  stress  components  must  be 
equal.  The  same  would  be  true  if  the  fluid  were  viscous. 

However,  the  classical  approach  with  an  inviscid  fluid-solid  interface  is  to  employ 
perfect-slip  boundary  conditions  that  allows  discontinuity  in  planar  displacement  com¬ 
ponents.  That  is,  the  radial  component  of  displacement  of  the  fluid  and  solid  must  be 
equal  at  the  interfaces,  however  the  circumferential  and  longitudinal  components  are 
discontinuous.  The  three  surface  stresses  must  also  be  equal.  However,  we  earlier  ob¬ 
served  that  the  stress  tensor  for  an  inviscid  fluid  is  isotropic.  Therefore,  the  shear  stress¬ 
es  in  the  inner  and  outer  fluids  are  identically  zero. 
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These  boundary  conditions  can  be  written  as 


[V  ''re*  r=  a  =  [“Pi’  =  a»  (2.61) 

and 

[V.W“,l,.b=  [-P2.0.0.>‘D,.b-  (2.62) 

Equations  2.61  and  2.62  represent  a  system  of  eight  linear,  algebraic  equations  in 
unknown  wave  amplitudes  that  in  matrix  form  are 

[L]  {c}  =  {0}  .  (2.63) 

where  [L]  is  a  8  x  8  matrix  whose  components  are  found  from  equations  2.35  and  2.39 
to  2.41,  and  from  equations  2.57  to  2.60,  and  where  {c}  is  an  8  x  1  column  vector  of 
the  unknown  amplitude  coefficients  Aj,  Bj,  A2,  B2,  A3,  B3,  D^,  and  D2.  The  compo¬ 
nents  of  [L]  are  defined  in  appendix  A. 

The  solution  of  equation  2.63  is  nontrivial  when  the  determinant  of  the  coefficients 
of  the  wave  amplitudes  {c}  vanishes,  i.e.  when, 

|L|  =0.  (2.64) 
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Equation  2.64  is  the  frequency  equation  of  the  coupled  system  consisting  of  a  trans¬ 
versely  isotropic  cylinder  filled  with  and  immersed  in  an  inviscid  fluid.  To  the  author  s 
knowledge,  this  frequency  equation  has  not  been  developed  before  within  the  frame¬ 
work  of  the  theory  of  elasticity. 
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3.  COMPUTATIONAL  RESULTS  FOR  ISOTROPIC  CYLINDERS 


We  begin  our  investigations  by  specializing  the  coupled  system  for  the  case  of  an 
isotropic  cylinder  filled  with  and  immersed  in  an  inviscid  fluid,  and  the  equations  of 
motion  of  the  n=l,  nonaxisymmetric  modes  of  wave  propagation.  The  cut-off  frequen¬ 
cies  are  computed  for  a  ratio  of  inner  to  outer  radius  s  ranging  from  0  to  0.9,  which  cor¬ 
responds  to  cylinders  ranging  from  solid  to  thin-walled.  The  frequency  spectra  are 
computed  for  relatively  thick  cylinders,  with  a  ratio  of  inner  to  outer  radius  of  0.8,  con¬ 
sisting  of  two  different  materials,  namely  a  soft  (linear)  rubber  and  steel.  The  relevant 
properties  of  the  two  solid  materials  and  of  the  inner  and  outer  fluids  are  provided  in 
appendix  C. 

Numerical  results  are  obtained  for  a  hollow  cylinder  in  a  vacuum,  a  hypothetical 
fluid  column,  a  fluid-filled  cylinder,  and  a  cylinder  that  is  fluid-filled  and  immersed  in 
fluid.  These  individual  cases  provide  a  logical  progression  to  systems  of  increasing  dif¬ 
ficulty. 

To  assist  in  the  interpretation  of  the  results  for  the  relatively  complex  case  of  a  cyl¬ 
inder  that  is  fluid-filled  and  immersed  in  fluid,  in  terms  of  the  results  for  simpler  con¬ 
figurations,  the  names  of  the  modes  are  consistent  from  configuration  to  configuration 
whenever  the  cut-off  frequencies  of  the  modes  can  be  approximately  traced  back  to  the 
corresponding  cut-off  frequencies  for  a  simpler  configuration.  The  results  for  the  iso¬ 
tropic  cylinders  will  be  used  as  an  aid  to  the  interpretation  of  the  results  for  a  cylinder 
of  a  transversely  isotropic  material. 
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The  numerical  results  are  presented  in  and  discussed  as  functions  of  normalized  fre¬ 
quency  Q ,  normalized  wavenumber  6 ,  normalized  phase  velocity  c ,  and  the  ratio  of 
the  cylinder’s  inner  to  outer  radius  s,  where  these  normalized  variables  are  defined  as 


a  =  inner  radius, 
b  =  outer  radius, 

s  =  ratio  of  inner  and  outer  radius, 
b 

h  =  b  -  a,  cylinder  thickness, 

Q  =  —,  normalized  frequency, 

5  =  kh,  normalized  wavenumber, 
c 

c  =  -^  =  —,  normalized  phase  velocity, 
0 


and  where 


is  the  shear  wave  velocity  in  the  material. 

3.1  ISOTROPIC  BEHAVIOR 

We  specialize  the  cylinder’s  constitutive  behavior  from  transversely  isotropic  to 
isotropic  by  reducing  the  number  of  elastic  constants  from  five  to  two  using  the  follow¬ 
ing  relationships 
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(3.2) 


C33  =  Cjj 
Ci3  =  Cj2 

C44  -  ^66  ~  o 


These  elastic  constants  are  written  in  tenns  of  engineering  constants  (also  known  as 
technical  constants)  as  follows, 


_  E(l-v) 

”  (l+v)(l-2v) 

Ev 

^12  (l+v)(l-2v) 

E 

^44  2(1+V)’ 


(3.3) 


where  E  is  the  Young’s  modulus  and  v  is  the  Poisson’s  ratio  of  the  isotropic  material. 
By  normalizing  c^  and  C12  to  C44,  these  elastic  constants  simplify  to 


2(l-v) 

"'ll  ■  (1-2v)''44 
2v 

-  (1-2v)''44. 


(3.4) 


Substitution  of  equations  3.2  and  3.4  into  equations  2.23,  2.26  and  2.38  results  in 
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■'44 


^11 


2  0(0  ,2 

q  -^- - k  = 

C44 


ik^  _k^ 

^  C44 

ik'n2  -k^- 


( 1  -  2v)  pco^ 
2(l-v)  C44 


(l-2v)co^  .2 
2(l-v) 


(3.5) 


3.2  A  HOLLOW  CYLINDER  IN  A  VACUUM 

We  first  consider  the  case  of  a  hollow  cylinder  with  no  external  or  internal  fluid 
loading.  For  this  set  of  boundary  conditions,  the  frequency  equation  is 


Lii 

^12 

Li3 

^14 

^15 

^16 

Lji 

L22 

^23 

^24 

^25 

^26 

L31 

^2 

^3 

L34 

^5 

^6 

^l 

^2 

^3 

L54 

L55 

L56 

L61 

^62 

^63 

L64 

L65 

^66 

^2 

^3 

L74 

^75 

^76 
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where  the  simplified  variables  in  equations  3.5  are  substituted  appropriately  into  the 
components  of  the  determinant. 

3.2.1  Cut-off  frequencies 

We  begin  by  solving  for  the  cut-off  frequencies  of  equation  3.6.  When  the  wave- 
number  is  zero,  the  variables  in  equations  3.5  further  reduce  to 


yl  CO 
Cnr* 


yl  (1-2v)q) 


2  CD 

q 


CO 

ikTlj  -^-2 
Cnn 


ikTi2  0 


(3.7) 


and  L23»  L24,  L25,  L26,  Lg3,  L54,  L55,  and  are  zero.  Equation  3.6  separates  into  the 
prcxiuct  of  two  sub-determinants 


DetjDet2  =  0 


(3.8) 


where 
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Detj  = 


Li3  Li4  Li5  Ljg 
^33  ^34  ^35  ^36 
^53  ^54  ^55  ^56 

L73  L74  L75 


(3.9) 


and 


Det2 


^21  ^22 
^62 


(3.10) 


Equation  3.8  is  satisfied  if  either  Detj  or  Det2  is  equal  to  zero.  The  case  Detj  =  0 
corresponds  to  plane  strain  vibrations  (Gazis  1958)  in  which  the  displacement  compo¬ 
nent  Uz  is  equal  to  zero.  The  case  Det2  =  0  corresponds  to  axial  (or  longitudinal)  shear 
vibration  involving  only  longitudinal  displacements  u^.  When  the  variables  in  equa¬ 
tions  3.7  are  substituted  into  equation  3.10,  we  find  that  the  only  material  variable  in 
Det2  is  the  shear  wave  velocity  in  the  material.  Therefore,  the  normalized  cut-off  fre¬ 
quencies  for  the  axial  shear  modes  are  the  same  for  any  isotropic  material.  Although  the 
plane  strain  and  axial  shear  modes  are  uncoupled  in  the  limiting  case  of  zero  wavenum¬ 
ber  (corresponding  to  rigid  body  motion),  these  modes  are  coupled  at  all  other  wave- 
numbers. 
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There  is  an  infinite  number  of  cut-off  frequencies  that  are  the  solutions  of  setting 
Det  j  or  Det2  to  zero.  For  the  case  of  a  hollow  cylinder  in  a  vacuum,  we  have  solved  for 
the  first  four  cut-off  frequencies  for  the  plane  strain  modes  and  the  first  two  cut-off  fre¬ 
quencies  for  the  axial  shear  modes.  It  should  be  noted  that  the  cut-off  frequency  for  the 
first  plane  strain  mode  is  zero.  A  plot  of  these  cut-off  frequencies  for  a  Poisson’s  ratio 
of  0.3  (representative  of  steel)  and  for  a  Poisson’s  ratio  of  0.45  (representative  of  rub¬ 
ber),  as  functions  of  the  ratio  of  the  inner  radius  and  outer  radius  of  the  cylinder,  is 
shown  in  figure  1.  (Note  that  figures  1-10  are  located  at  the  end  of  chapter  3,  starting 
on  page  65). 

The  cut-off  frequencies  are  normalized  to  the  shear  wave  velocity  in  each  material. 
Using  the  material  properties  for  rubber  and  for  steel,  provided  in  appendix  C,  along 
with  equation  3.1,  the  shear  wave  velocities  in  the  two  materials  are 

c  j  =  72  m/sec,  for  rubber 
Cy  =  3218  m/sec,  for  steel. 


The  normalized  cut-off  frequencies  for  the  axial  modes  are  the  same  for  either  ma¬ 
terial  type,  as  was  discussed  above.  The  normalized  cut-off  frequencies  for  the  second 
and  third  plane  strain  modes  are  similar  for  both  materials,  with  those  for  rubber  being 
somewhat  greater  than  those  for  steel.  The  normalized  cut-off  frequencies  for  the  sec¬ 
ond  and  third  plane  strain  modes,  of  the  two  materials,  become  asymptotically  close  as 
the  thickness  of  the  cylinder  wall  decreases.  The  normalized  cut-off  frequencies  for  the 
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fourth  plane  strain  mode  are  quite  different  for  the  two  materials,  particularly  as  the 
wall  thickness  increases.  As  the  wall  thickness  decreases,  the  normalized  cut-off  fre¬ 
quencies  of  the  fourth  plane  strain  modes,  for  both  materials,  become  essentially  con¬ 
stant,  but  not  equal. 

The  results,  shown  in  figure  1,  for  v  =  0.3 ,  agree  with  those  of  Armenakas  et  al. 
(1969),  Gazis  (1959),  and  Kumar  and  Stephens  (1972).  Kumar  and  Stephens  (1972) 
also  investigated  the  variation  of  cut-off  frequencies  as  functions  of  Poisson  s  ratio  as 
Poisson’s  ratio  was  varied  from  0  to  0.5,  for  the  case  of  a  ratio  of  inner  radius  to  outer 
radius  of  0.5.  Their  results  showed  that,  in  general,  the  cut-off  frequencies  increased 
with  an  increasing  Poisson’s  ratio. 

The  cut-off  frequencies  for  s=0.8  will  be  used  to  identify  the  modes  in  the  frequency 
spectrum  for  each  material. 

3.2.2  Frequency  spectra 

Next  we  develop  the  branches  of  the  frequency  spectrum  for  each  material  by  nu¬ 
merically  searching  for  the  frequencies,  which  are  solutions  of  equation  3.6,  for  real 
wavenumbers  from  0  to  k.  Figures  2  and  3  are  plots  of  the  frequency  spectrum  for  the 
three  lowest  modes  (first  two  plane  strain  modes  and  first  axial  mode)  in  a  cylinder 
made  of  a  soft  rubber  and  steel,  respectively. 

In  both  figures,  the  first  branch,  labeled  Sj,  corresponds  to  the  first  plane  strain 
mode.  This  mode  is  a  flexural  mode  with  large  radial  displacements  and  a  zero  cut-off 
fi^uency.  The  second  branch,  labeled  S2,  corresponds  to  an  axial  shear  mode  with 
large  axial  displacements  and  whose  cut-off  frequency  is  the  lowest  solution  of  setting 
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Det2  equal  to  zero.  For  both  materials  the  normalized  cut-off  frequency  for  S2  is 


Q  =  0.223  . 


The  third  branch,  labeled  S3,  corresponds  to  the  second  plane  strain  mode.  This 
mode  is  a  “breathing  mode”  with  large  tangential  displacements.  The  cut-off  frequency 
for  S3  is  the  lowest,  non-zero  solution  of  setting  Det^  equal  to  zero.  The  normalized  cut¬ 
off  frequencies  for  the  two  materials  are 

Q  s  0.5953,  for  rubber 
Q  =  0.5274,  for  steel. 

For  both  materials,  S2  and  S3  have  group  velocities  that  are  vanishingly  small  in  the 
neighborhood  of  the  cut-off  frequencies. 

The  low  wavenumber  behavior  of  the  Sj,  S2,  and  S3  modes  corresponds  to  that  con¬ 
tained  in  a  bending  shell  theory  in  which  the  radial  component  of  displacement  is  as¬ 
sumed  constant  through  the  thickness,  while  the  axial  and  tangential  components  are 
assumed  to  vary  linearly  through  the  thickness  (Armenakas,  1969). 

By  comparing  figures  2  and  3,  we  can  examine  the  differences  in  the  wave  propa¬ 
gation  characteristics  of  hollow  cylinders  of  rubber  and  of  steel.  For  instance,  the  nor¬ 
malized  phase  velocity  of  S 1,  for  both  rubber  and  steel,  is  essentially  identical  in  the  low 
wavenumber  domain.  At  higher  wavenumbers,  the  normahzed  phase  velocity  of  the 
flexural  wave  of  a  steel  cylinder  is  less  than  that  of  a  rubber  cylinder  of  comparable 
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thickness.  The  normalized  cut-off  frequency  and  normalized  phase  velocity,  over  all 
wavenumbers,  of  S2  is  identical  for  both  materials  and  is  asymptotic  to  the  shear  wave 
velocity  in  each  material  at  high  wavenumbers.  The  normalized  cut-off  frequency  and 
normalized  phase  velocity  of  S3  for  the  steel  cyhnder  is  lower  than  for  the  rubber  cyl¬ 
inder,  over  all  wavenumbers  shown. 

The  frequency  spectrum  for  the  steel  cylinder  compares  very  well  with  that  of  pre¬ 
vious  investigators  (Gazis  1959,  Armenakas  et  al.  1969,  and  Kumar  and  Stephens 
1972).  These  numerical  results  are  included  here  for  completeness  and  as  an  aid  to  the 
interpretation  of  the  results  for  more  complicated  geometries.  The  frequency  spectrum 
for  the  thick  rubber  cylinder,  to  the  author’s  knowledge,  has  not  been  presented  before. 
33  A  HYPOTHETICAL  FLUID  COLUMN 

Earlier  studies  by  Kumar  (1971, 1972),  Chandra  and  Kumar  (1976),  and  Sinha  et 
al.  (1992)  have  shown  that  the  presence  of  an  inner  fluid  may  introduce  modes  of  wave 
propagation  that  are  not  present  in  a  hollow  cylinder  in  a  vacuum.  These  additional 
modes  originate  from  modes  in  the  fluid  column.  Therefore,  we  next  consider  the  cut¬ 
off  frequencies  and  frequency  spectrum  of  a  hypothetical,  cylindrical  fluid  column. 

This  information  will  be  useful  to  the  interpretation  of  the  frequency  spectrum  of 
the  case  of  a  fluid-filled  cylinder. 

3  J.l  Pressure-release  boundary  conditions 

The  case  of  a  cylindrical  fluid  column  with  a  pressure-release  boundary  condition 
requires  the  pressure  at  the  inner  fluid-solid  interface  to  vanish  while  the  radial  dis¬ 
placement  remains  arbitrary.  Since  the  pressure  must  be  zero  for  all  0,  z,  and  t,  equation 
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2.58,  at  the  inner  radius  and  n=l,  reduces  to 


Wi(44a)=0.  (3.11) 

or 

^^a  =  XpX2,  X3, ...,  (3.12) 

where  Xj,  X2,  X3, ...  are  the  zeroes  of  Wj  .  If  Wj  (^4 a)  =  Ij  (^4a)  then 
Xp  X2,  X3, . . . ,  are  all  complex  (Abramowitz  and  Stegun  1972)  and,  from  equation  3.12, 
( 44a)  is  complex.  The  inner  radius  a  is  always  a  real,  positive  quantity,  therefore  com¬ 
plex  (^4a)  implies  complex  ^4 .  However,  as  was  discussed  in  chapter  2,  if  ^4  is  com¬ 
plex  then  Wj  (^4a)  =  Jj  (^4a)  .  This  is  contradictory.  So  Wj  (^4a)  ^  .  If 

Wj  (^4a)  =  (^4a)  then  Xp  X2,  X3, ... ,  are  all  real.  This  means  (^4a)  is  real  and, 

since  the  inner  radius  a  is  always  a  real,  positive  quantity,  ^4  is  real  and  is  equal  to  a . 
Therefore,  from  equation  3.11,  we  have 

Ji(aa)=0,  (3.13) 

and,  from  equation  3.12,  we  have 

aa  =  3.832, 7.016, 10.173, ....  (3.14) 
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Equation  3. 14  is  the  frequency  equation  of  a  hypothetical  fluid  column  with  a  pressure- 
release  boundary  condition. 

3.3.2  Displacement-free  boundary  conditions 

The  case  of  a  cylindrical  fluid  column  with  radial  displacement-free  boundary  con¬ 
ditions  requires  that  u^^  must  vanish  at  the  inner  radius  of  the  cylinder  while  the  pressure 
remains  arbitrary.  This  must  be  tme  for  all  0,  z,  and  r,  and  equation  2.57 ,  at  the  inner 
radius  and  n=l,  reduces  to 


[W,(44r)  1  _  =  0.  (3.15) 

Using  a  similar  argument  as  above,  Wj  (^4a)  *  Ij  (^43)  because  the  zeros  of 
Ij  (^4r)  ^  are  all  complex  (Abramowitz  and  Stegun  1972).  If  Wj  (^4 a)  =  J j  (^4a) 
then  the  roots  of  [  J,  (£.r)  ]  are  all  real,  =  a ,  and  equation  3.15  becomes 

Jj  (aa) -aaJ2(aa)  =0,  (3.16) 

or, 

aa  =  1.841, 5.331,  8.536, ... .  (3.17) 

Equation  3.17  is  the  frequency  equation  of  a  hypothetical  fluid  column  with  a  displace¬ 
ment-free  boundary  condition. 

Equations  3.14  and  3.17  represent  families  of  hyperbolas  with  finite  cut-off  fre- 


53 


quencies.  The  first  two  cut-off  frequencies,  for  each  of  the  two  boundary  conditions  and 
for  s=0.8,  normalized  to  the  shear  wave  velocity  in  rubber,  are 

Q  =  19.98,  36.58,  pressure-release  boundary  conditions 
Q  =  9.61,  27.8,  displacement-free  boundary  condtions. 

The  first  two  cut-off  frequencies,  for  each  of  the  two  boundary  conditions  and  for  s=0.8, 
normalized  to  the  shear  wave  velocity  in  steel,  are 

Q  =  0.9465, 0.8176,  pressure-release  boundary  conditions 
=  0.2146,  0.6213,  displacement-free  boundary  conditions. 

3.3.3  Frequency  spectra 

The  first  two  branches,  corresponding  to  each  of  the  two  sets  of  boundary  condi¬ 
tions  and  normalized  to  the  shear  wave  velocity  in  rubber  and  of  steel,  are  shown  in  fig¬ 
ures  4  and  5,  respectively.  The  normalization  to  the  shear  wave  velocity  in  each  of  the 
materials  facilitates  the  comparison  of  these  results  to  other  results.  The  curves  labeled 
l^and  I2  represent  solutions  of  equations  3.14,  while  the  curves  labeled  and  L2  rep¬ 
resent  solutions  of  equations  3.17. 

From  a  comparison  of  figures  2  and  4,  we  note  that  the  fluid  modes,  figme  4,  exist 
at  wavenumbers  and  frequencies  greater  than  the  solid  modes  of  a  rubber  cylinder,  fig¬ 
ure  2.  This  means  that  the  fluid  modes  will  not  couple  with  the  solid  modes,  at  least  in 
the  wavenumber  and  frequency  domains  of  interest.  However,  from  a  comparison  of 
figures  3  and  5,  we  note  that  the  fluid  modes,  figure  5,  exist  within  the  same  wavenum- 
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ber  and  frequency  domains  as  the  solid  modes  of  a  steel  cylinder,  figure  3.  The  impli¬ 
cations  of  this  overlap  or  intersection  of  fluid  and  solid  modes  will  be  discussed  further 
in  the  next  section. 

Sinha  et  al.  (1992)  investigated  the  dispersion  curves  for  the  axisymmetric  wave 
propagation  in  a  hypothetical  fluid  column,  with  stress-free  and  with  displacement-free 
boundary  conditions,  and  normalized  their  results  to  the  shear  wave  velocity  in  steel. 
They  found  that  the  lowest  order  dispersion  curve  for  the  displacement-free  boundary 
condition  had  a  zero  cut-off  frequency.  Kumar  (1971)  derived  the  frequency  equation 
of  the  n=l  modes  of  nonaxisymmetric  wave  propagation  in  a  hypothetical  fluid  column, 
for  stress-free  and  for  displacement-free  boundary  conditions.  However,  he  did  not  de¬ 
velop  the  dispersion  curves. 

3.4  A  FLUID-FILLED  CYLINDER 

The  next  case  to  be  investigated  is  that  of  a  fluid-filled  cylinder.  For  this  set  of 
boundary  conditions,  the  frequency  equation  is 
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In  the  limiting  case  of  zero  wavenumber,  this  frequency  equation  separates  into  the 
product  of  two  sub-determinants 


DetjDet2  =  0 


(3.19) 


where 
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and  where 


Det2 


^21  ^2 
^62 


(3.21) 


is  the  same  as  equation  3.10. 

In  the  case  of  a  hollow  cylinder  without  any  fluid  loading,  we  noted  that  the  only 
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material  variable  in  Det2  was  the  shear  wave  velocity  in  the  solid  material.  Inspection 
of  the  components  of  equation  3.21  shows  that  Det2  is  also  independent  of  the  internal 
fluid  parameters.  This  means  that  the  normalized  cut-off  frequencies  for  the  axial  shear 
modes,  found  from  setting  3.10  or  3.21  equal  to  zero,  are  the  same  for  empty  and  for 
fluid-filled  cylinders  of  any  isotropic  material. 

3.4.1  Cut-off  frequencies  for  the  rubber  cylinder 

The  normalized  cut-off  frequencies  for  the  Si,  S2,  and  S3  modes  for  the  fluid-filled 
rubber  cylinder  are 

flsO,  0.223,  0.431, 

respectively. 

As  expected,  the  presence  of  the  internal  fluid  has  not  affected  the  normalized  cut¬ 
off  frequency  of  the  S2  mode.  However,  the  presence  of  the  internal  fluid  has  reduced 
the  cut-off  frequency  of  the  S3  mode  by  approximately  28%.  No  additional  cut-off  fre¬ 
quencies,  which  can  be  traced  back  to  fluid  modes,  exist  within  the  normalized  frequen¬ 
cy  band  of  0  to  it.  The  absence  of  any  additional  cut-off  frequencies  means  that  the 
internal  fluid  does  not  introduce  additional  modes  in  the  coupled  system,  consisting  of 
the  inner  fluid  and  the  cylinder,  within  the  normalized  frequency  band  of  interest. 

3.4.2  Frequency  spectrum  for  the  rubber  cylinder 

The  frequency  spectrum  for  the  fluid-filled  rubber  cylinder  is  shown  in  figure  6. 
From  a  comparison  of  figure  6  to  figure  2,  we  note  that  the  phase  velocity  of  the  S^ 
mode  of  the  fluid-filled  cylinder  is  less  than  that  of  the  Si  mode  of  the  empty  cylinder 
for  all  wavenumbers.  For  instance,  at  5  =  0.271,  the  phase  velocity  of  the  Si  mode  is  re- 
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duced  approximately  36%  when  fluid  is  added  to  the  inside  of  the  cylinder. 

The  phase  velocity  of  the  S2  mode  is  reduced  somewhat  for  d/n  <  0.3 .  At  higher 
wavenumbers,  the  phase  velocity  of  this  mode  is  unaffected  by  the  fluid  and  is  asymp¬ 
totic  to  the  shear  wave  velocity  in  mbber.  The  phase  velocity  of  the  S3  mode  is  reduced 
over  all  wavenumbers  as  compared  to  the  phase  velocity  of  the  hollow  cylinder  without 
an  internal  fluid. 

3.4.3  Cut-off  frequencies  for  the  steel  cylinder 

The  normalized  cut-off  frequencies  for  the  Sj,  S2,  Lj,  S3,  and  L2  modes  for  a  fluid- 
filled  steel  cylinder  are 

Q  s  0, 0.223, 0.228, 0.516, 0.644, 

respectively. 

The  normalized  cut-off  frequencies  of  each  of  these  modes  can  be  traced  back  to  the 
normalized  cut-off  frequencies  of  the  modes  for  each  of  the  separate  components  of  the 
coupled  system.  For  instance,  the  cut-off  frequency  for  the  S3  mode  is  somewhat  re¬ 
duced  from  the  cut-off  frequency  for  the  “breathing”  mode  of  the  empty  cylinder.  The 
cut-off  frequencies  for  the  modes  labeled  Lj  and  L2  can  be  similarly  traced  back  to  the 
cut-off  frequencies  of  the  hypothetical  fluid  column  with  the  radial  displacement-free 
boundary  condition.  As  expected,  the  cut-off  frequency  of  the  S2  mode  has  been  unaf¬ 
fected  by  the  presence  of  the  fluid. 

The  presence  of  the  additional  cut-off  frequencies  for  the  L^and  L2  modes  means 
that  the  internal  fluid  has  introduced  two  more  modes  in  the  coupled  system,  consisting 
of  a  the  inner  fluid  and  the  cylinder,  at  least  within  the  normalized  frequency  band  of  0 
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t0  7t. 


3.4.4  Frequency  spectrum  for  the  steel  cylinder 

The  frequency  spectrum  for  the  fluid-filled  steel  cylinder  is  shown  in  figure  7. 
These  results  are  dramatically  different  from  those  shown  in  figure  6.  The  fluid-based 
modes  and  the  solid-based  modes  are  greatly  perturbed  in  the  coupled  system  as  evi¬ 
denced  by  large  changes  in  the  phase  velocities  from  those  of  the  individual  compo¬ 
nents  of  the  system  and  by  an  alternating  “weak  coupling”  between  the  S3  and  modes 

and  between  the  S3  and  modes.  Weak  coupling  of  modes  is  evidenced  by  regions  in 

the  frequency  spectrum  where  branches  come  close  together  (Armenakas  et  al.  1969). 

The  phenomenon  of  coupling  between  fluid  and  solid  modes  was  observed  for  the 
cases  of  axisymmetric  motion  of  a  fluid-filled  steel  cylinder  by  Sinha  et  al.  (1992)  and 
for  the  case  of  nonaxisymmetric  motion  of  a  fluid-filled  steel  cylinder  by  Kumar 
(1971). 

A  comparison  of  the  frequency  spectrum  of  a  fluid-filled  steel  and  hard  rubber  cyl¬ 
inders  was  made  by  Fuller  and  Fahy  (1982)  on  the  basis  of  a  thin-shell  theory.  In  their 
analysis  the  ratio  of  the  cylinder  wall  thickness  to  the  mean  radius  was  0.05,  and  the 
shear  wave  velocity  in  the  hard  rubber  was  approximately  7 85  m/sec.  The  results  shown 
in  figures  6  and  7  are  for  a  ratio  of  the  cylinder  wall  thickness  to  the  mean  radius  of 
approximately  0.22,  and  for  a  shear  wave  speed  in  the  soft  rubber  of  approximately  72 
m/sec. 

To  the  author’ s  knowledge,  this  is  the  first  time  that  a  direct  comparison  of  the  wave 
propagation  characteristics  of  a  thick,  fluid-filled  steel  and  soft  robber  cylinder  has  been 
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made  on  the  basis  of  the  theory  of  elasticity. 

3.5  A  CYLINDER  THAT  IS  FLUID-FILLED  AND  IMMERSED  IN  FLUID 

The  final  case  to  be  investigated,  in  this  chapter,  is  that  of  a  fluid-filled  cylinder  that 
is  immersed  in  an  infinite  inviscid  fluid.  For  this  set  of  boundary  conditions,  the  fre¬ 
quency  equation  is 
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In  the  limiting  case  of  zero  wavenumber,  this  frequency  equation  separates  into  the 
product  of  two  sub-determinants 


DetjDet2  =  0 


(3.23) 


where 
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and  where 


Det2 


^21  ^22 
^62 


(3.25) 


is  the  same  as  equations  3.10  and  3.21.  Inspection  of  the  components  of  equation  3.25 
shows  that  Det2  is  also  independent  of  the  external  fluid  parameters.  Therefore,  the  nor¬ 
malized  cut-off  frequencies  for  the  axial  shear  modes,  found  from  setting  equations 
3.10, 3.21,  or  3.25  equal  to  zero,  are  the  same  for  any  isotropic  cylinders  that  are  empty, 
fluid-filled,  or  fluid-filled  and  immersed  in  fluid. 

3.5.1  Cut-off  frequencies  for  the  rubber  cylinder 

The  normalized  cut-off  frequencies  for  the  Sj,  S2,  and  S3,  modes  for  a  rubber  cyl¬ 
inder  that  is  fluid-filled  and  immersed  in  fluid  are 

12  =  0,0.223, 0.431, 
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respectively. 

The  cut-off  frequencies  for  the  S3  mode  is  essentially  unaffected  by  the  addition  of 
the  external  fluid  loading.  All  the  cut-off  frequencies  are  real. 

3  J.2  Frequency  spectrum  for  the  rubber  cylinder 

The  frequency  spectrum  for  the  rubber  cyhnder  that  is  fluid-filled  and  immersed  in 
fluid  is  shown  in  figure  8.  From  a  comparison  of  figures  6  and  8,  we  note  that  the  low 
wavenumber  behavior  of  the  Sj  and  S3  modes  is  virtually  unaffected  by  the  addition  of 
the  external  fluid  loading.  There  is  some  reduction  in  phase  velocity  of  the  S  j  mode  at 
the  higher  wavenumbers.  The  reduction  in  phase  velocity  of  the  S3  mode  at  the  higher 
wavenumbers  is  more  pronounced.  The  S2  mode  is  not  affected  by  the  external  fluid 
loading.  It  should  be  noted,  that  the  branches  are  real.  This  means  that  the  external  fluid 
has  not  introduced  any  damping  into  the  modes  of  wave  propagation  in  the  coupled  sys¬ 
tem. 

3^.3  Cut-off  frequencies  for  the  steel  cylinder 

The  normalized  cut-off  frequencies  for  the  Sj,  S2,  L^,  S3,  and  L2  modes  for  a  steel 
cylinder  that  is  fluid-filled  and  immersed  in  fluid  are 

Q  =  0, 0.223, 0.228  -  0.0014i,  0.515  -  O.OlOi,  0.6433-0.00351, 
respectively. 

The  most  notable  effect  of  the  external  fluid  loading  on  the  cut-off  frequencies  of 
the  coupled  system  is  that  the  cut-off  frequencies  for  the  L^,  S3,  and  L2  modes  are  com¬ 
plex.  At  zero  wavenumber,  these  modes  are  rigid  body  modes.  Therefore,  the  presence 
of  the  external  fluid  has  introduced  damping  in  these  modes  of  vibration. 
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3.5.4  Frequency  spectrum  for  the  steel  cylinder 

In  figure  9  we  show  the  real  part  of  the  frequency  spectrum  for  the  steel  cylinder 
that  is  fluid-filled  and  immersed  in  fluid.  From  a  comparison  of  figures  9  and  7  we  find 
that  the  real  part  of  branches  are  very  similar  to  the  branches  for  the  fluid-filled  steel 
cylinder.  That  is,  except  some  very  small  changes  in  the  phase  velocity  of  the  Sj  mode, 
the  presence  of  the  external  fluid  has  had  little  effect  on  the  real  part  of  the  cut-off  fre¬ 
quencies  or  the  phase  velocities  of  any  of  the  modes. 

In  figure  10  we  show  the  imaginary  part  of  the  frequency  spectrum.  The  presence 
of  the  imaginary  part  of  a  complex  frequency  indicates  that  the  external  fluid  has  intro¬ 
duced  damping  in  the  coupled  system.  What  is  most  important  to  note  in  figure  10  is 
the  nontrivial  behavior  of  the  imaginary  term  as  a  function  of  wavenumber  and  of  mode 
type.  At  zero  wavenumber,  the  Sj  and  S2  modes  are  undamped  and  the  cut-off  frequen¬ 
cies  are  unaffected  by  the  presence  of  the  internal  or  external  fluids.  The  Sj  mode  re¬ 
mains  lightly  damped  or  undamped  over  the  wavenumber  domain  of  interest.  The  S3 
mode  has  the  highest  amount  of  damping  at  zero  wavenumber  while  the  damping  in  the 
L2  mode  is  somewhat  greater  than  the  damping  in  the  mode.  This  means,  as  rigid 
body  modes,  only  the  Sj  (harmonic  transverse  vibration)  and  S2  modes  (harmonic  lon¬ 
gitudinal  vibration)  will  exist  for  all  time. 

Attenuation  of  displacements,  particle  velocities,  stresses,  etc.,  due  to  the  damping 
in  the  system,  can  be  examined  by  rewriting  the  complex  exponential  as 


where  x  is  the  real  part  and  y  is  the  imaginary  part  of  the  complex  frequency.  If  y  is 
positive,  the  cylinder  displacements  (for  instance)  increase  without  bounds  with  in¬ 
creasing  time.  This  is  physically  unrealistic.  Therefore,  the  imaginary  term  must  be 
negative,  which  is  true  for  our  results. 

Chandra  and  Kumar  (1976)  investigated  the  axisymmetric  modes  of  wave  propaga¬ 
tion  in  a  brass  cylinder  that  was  fluid-filled  and  immersed  in  fluid.  They  found  that  the 
percent  reduction  in  phase  velocity  of  the  various  modes  increased  with  a  decrease  in 
the  wall  thickness  of  the  cylinder.  Sinha,  et  al.  (1992)  investigated  the  axisymmetric 
modes  of  wave  propagation  in  a  steel  cylinder  that  was  fluid-filled  and  immersed  in  flu¬ 
id.  They  found  the  presence  of  the  external  fluid  introduced  damping  in  the  system 
through  the  presence  of  complex  ftequencies  as  solutions  of  the  frequency  equation. 

To  the  author’s  knowledge,  this  dissertation  represents  the  first  time  the  frequency 
spectrum  of  the  n=l  modes  of  nonaxisymmetric  wave  propagation  in  isotropic  cylin¬ 
ders  (of  any  material),  with  internal  and  external  fluid  loading,  has  been  studied  within 
the  framework  of  the  theory  of  elasticity. 


64 


cn 

3 


52 

o 


0) 

c 

c 


o 

_o 

(3 

DC 


Figure  1.  Comparison  of  the  normalized  cut-off  frequencies  of  the  steel  cylinder  and  of 
the  rubber  cylinder,  as  functions  of  the  ratio  of  the  inner  to  the  outer  radius 
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Figure  2.  Frequency  spectrum  of  the  rubber  cylinder  in  a  vacuum 
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Figure  3.  Frequency  spectrum  of  the  steel  cylinder  in  a  vacuum 
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Figure  4.  Frequency  spectrum  of  the  hypothetical  fluid  column  normalized  to  the  shear 
wave  velocity  in  rubber 
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Figure  5.  Frequency  spectrum  of  the  hypothetical  fluid  column  normalized  to  the  shear 
wave  velocity  in  steel 
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Figure  6.  Frequency  spectrum  of  the  fluid-filled  rubber  cylinder 
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Figure  7.  Frequency  spectrum  of  the  fluid-filled  steel  cylinder 
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Figure  8.  Frequency  spectrum  of  the  rubber  cylinder  that  is  fluid-filled  and  immersed 
in  fluid 
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Figure  9.  Real  part  of  the  frequency  spectrum  of  the  steel  cylinder  that  is  fluid-filled  and 
immersed  in  fluid 
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Figure  10.  Imaginary  part  of  the  frequency  spectrum  of  the  steel  cylinder  that  is  fluid 
filled  and  immersed  in  fluid 
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4.  COMPUTATIONAL  RESULTS  FOR  A  HOMOGENEOUS, 
TRANSVERSELY  ISOTROPIC  CYLINDER 

In  this  chapter,  the  n=l  modes  of  nonaxisymmetric  wave  propagation  in  a  fluid- 
loaded,  homogeneous,  transversely  isotropic  cylinder  are  investigated.  The  cut-off  fre¬ 
quencies  are  computed  for  a  ratio  of  cylinder  inner  to  outer  radius  s  ranging  from  0  to 
0.9.  The  frequency  spectrum  is  computed  for  a  cylinder  with  a  ratio  of  inner  to  outer 
radius  of  0.8  and  consisting  of  a  fiber-reinforced  matrix  material  in  which  the  fibers  are 
aligned  parallel  to  the  longitudinal  axis  of  the  cylinder. 

The  matrix  material  is  a  soft  rubber  with  the  same  material  properties  as  in  the  case 
of  the  rubber  isotropic  cylinder.  The  fibers  are  nylon  cords  and  the  fiber  volume  fraction 
is  10%.  The  material  properties  of  the  rubber  matrix  material  and  nylon  cords  are  pro¬ 
vided  in  appendix  C. 

Numerical  results  for  the  same  sequence  of  cases  of  boundary  conditions  are  ob¬ 
tained  in  this  chapter  as  were  obtained  in  chapter  3.  These  cases  are  a  hollow  cylinder 
in  a  vacuum,  a  hypothetical  fluid  column  with  pressure-release  and  with  displacement- 
free  boundary  conditions,  a  fluid-filled  cylinder,  and  a  cylinder  that  is  fluid-filled  and 
immersed  in  fluid. 

The  numerical  results  are  presented  in  and  discussed  as  functions  of  normalized  fre¬ 
quency  Q ,  normalized  wavenumber  8 ,  normalized  phase  velocity  Q/ 5 ,  and  the  ratio 
of  the  cylinder’s  inner  to  outer  radius  s.  These  normalized  variables  were  defined  in 
chapter  3,  but  are  redefined  here  for  completeness  as 
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a  =  inner  radius, 
b  =  outer  radius, 

s  =  ratio  of  inner  and  outer  radius, 
b 

h  =  b  -  a,  cylinder  thickness, 

Q  =  —,  normalized  frequency, 

Cj 

5  =  kh,  normalized  wavenumber, 

Q  c 

c  =  =  — ,  normalized  phase  velocity, 

0  Cnf 


where  we  define  the  shear  wave  velocity  in  the  transversely  isotropic  material  as 


c 


2 

T 


P  ■ 


(4.1) 


4,1  EFFECTIVE  ELASTIC  CONSTANTS 

For  this  investigation,  we  use  the  definitions  derived  by  Jones  (1975)  for  the  effec¬ 
tive  engineering  constants  of  the  fiber-reinforced  cylinder.  These  definitions  are 
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where  E  ,G  ,v  ,p  ,  and  V„  are  the  Young’s  modulus,  shear  modulus,  Poisson’s 
ratio,  density  and  volume  fraction  of  the  matrix  material,  and  where 
Ep  Gp  Vp  Pp  and  are  the  Young’s  modulus,  shear  modulus,  Poisson’s  ratio,  densi¬ 
ty  and  volume  fraction  of  the  fiber,  respectively,  and  where 
B-j-  =  the  transverse  bulk  modulus, 

Gl  =  the  longitudinal  shear  modulus, 

Gx  =  the  transverse  shear  modulus, 

El  =  the  longitudinal  Young’s  modulus. 

Ex  =  the  transverse  Young’s  modulus, 

Vl  =  the  longitudinal  Poisson’s  ratio, 

Vx  =  the  transverse  Poisson’s  ratio. 

In  these  definitions,  “transverse”  refers  to  properties  in  the  plane  of  isotropy,  and  “lon¬ 
gitudinal”  refers  to  properties  along  the  material  preferred  direction. 

The  transverse  bulk  modulus  is  defined  by  Hashin  (1979)  as 
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(4.3) 


Substitution  of  the  material  properties  of  the  matrix  material  (soft  rabber)  and  the 
reinforcing  fibers  (nylon  cord),  along  with  a  10%  fiber  volume  fraction  of  fiber,  into 
equations  4.2  and  4.3,  results  in 

=  1.54  X  lo’  N/m^ 

Gl  =  5.75  X  10^  N/m^ 

=  5.75  X  10^  N/m^ 

El  =  5.64  X  10^  N/m^ 

E^,  =  1.67  X  10^  N/m^ 

Vl  =  0.42 
v.p  =  0.45 

p  =  1014  kg/m^. 

Sutcu  (1992)  analyzed  the  constraints  on  the  engineering  constants  for  a  transverse¬ 
ly  isotropic  material  and  reported  them  as 
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Gl>  Gi-j  E^j  Ej  >  0, 


-1  <Vy<  1, 


^  ( 1  —  Vj) 

L-  2  Ej^' 


(4.4) 


The  effective  engineering  constants  of  the  fiber-reinforced  cylinder,  computed  from  the 
definitions  by  Jones  (1975),  fall  within  the  allowable  bounds  of  equations  4.4. 

Substitution  of  the  effective  engineering  constants  into  the  definitions  of  the  elastic 
constants,  equations  2.11,  results  in 


c^i  =  2.12  X  10^ 

Ci2  =  9.69  X  10^ 

Cj3  =  1.3  X  10^ 

C33  =  5.74  X  10* 

C44  =  5.74  X  10^ 

C44  =  5.74  xl0^ 

When  these  elastic  constants  are  normalized  to  C44,  we  have 
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Cii-3-7  XC44 


c 

c 

c 

c 


12-  1-7  XC44 

13  -  ^-3  ^  ^44 
33  =  100  X 


66  -  ^44  • 


It  is  important  to  take  note  at  this  time  of  how  highly  anisotropic  the  fiber-rein¬ 
forced  composite  material  is  compared  to  naturally  occurring  hexagonal  crystal  mate¬ 
rials.  For  instance,  for  the  fiber-reinforced  material,  C33  >  27  x  Cjj ,  while,  from  Payton 
(1983),  we  find  that  the  elastic  constants  for  beryllium  are 


Cji  =  1-8  X  c^ 

^12  =  0.16  X 
Ci3  =  0.09  X  c^ 
^33  =  2-1  ^  C/]/) 

^66  —  ^  ^44  » 


and,  therefore,  C33  s  1.17  x  c^^ .  Similarly,  the  elastic  constants  for  ice  are 
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Cji  =  4.3xC44 
C  =  2. 1  X 

C13S1.6XC44 

C33  -  ^  ^^44 

^66—  ^^44  ’ 


and,  therefore,  C33  =  1.07  xc^i- 

It  is  also  interesting  to  compare  the  normalized  elastic  constants  for  the  two  isotro¬ 
pic  materials  used  in  chapter  3  to  those  for  the  fiber-reinforced  material.  Using  equation 
3.4  and  v  =  0.45  for  the  soft  rubber  cylinder  we  have 


^11  “  ^33  “  11  X  C44 
^12  “  <^13  “  ^^^44 
^66  “  ^44  • 

Similarly,  using  equation  3.3  and  v  =  0.3  for  the  steel  cylinder,  we  have 


Cji  —  ^33  —  3.5  X  c^^ 

^12  “  ^13  “  ^  ^44 

^66  “  ^44  • 

The  relationships  of  c^  and  C12  to  C44  in  the  fiber-reinforced  material  are  very  similar 
to  those  in  steel.  However,  the  relationships  of  C13  and  C33  to  C44  in  the  fiber-reinforced 
material  are  different  from  those  in  either  steel  or  rubber.  In  particular,  the  relationship 
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of  C33  to  C44  is  approximately  an  order  of  magnitude  greater  than  in  rubber.  Therefore, 
we  expect  the  numerical  results  for  the  transversely  isotropic  cylinder  to  be  different 
from  the  results  for  either  the  rubber  or  the  steel  cylinder. 

4.2  A  HOLLOW  CYLINDER  IN  A  VACUUM 

The  frequency  equation  for  this  set  of  boundary  equations  was  defined  by  equation 
3.6.  It  is  repeated  here  for  completeness  as 


|L| 


'11 

Li2 

^13 

Li4 

Li5 

^16 

'21 

^22 

L23 

^24 

L25 

^26 

'31 

^2 

L33 

^34 

^35 

^36 

'51 

^2 

^3 

L54 

L55 

L56 

'61 

^62 

^63 

L64 

^65 

L66 

'71 

^72 

L73 

L74 

^75 

^76 

(4.5) 


4.2.1  Cut-off  frequencies 

We  first  examine  the  solutions  of  equation  4.5  in  the  limiting  case  of  zero  wave- 
number.  When  the  wavenumber  is  zero,  equations  2.23, 2.26  and  2.38,  simplify  to 
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C  “ 
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^44  1(0 


■'44 


2  2 

J.2  ^ 

*2  ^  2 

^44  c^ 

<^66 

ikT]^  ->0 
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(0 


‘'ll  ^44 

II 

1 

3 

Q. 

‘^ll  “^^44 

^Ci3-t-C44^ 

IC44  1 

^Ci3  -t-C44^ 

2’ 


(4.6) 


and  L21,  L22>  L25,  L26’  ^61’  ^62.  ^65’  andLg^  are  zero.  Just  as  in  the  case  of  the  isotropic 
cylinders,  equation  4.6  separates  into  the  product  of  two  sub-determinants 


DetjDet2  =  0 


(4.7) 


where 


Detj  = 


^12  ^15 
^32  ^35 
^52  ^55 
^72  ^75 


(4.8) 
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and 


Detj 


L'23  ^24 
^63  ^64 


(4.9) 


Equation  4.7  is  satisfied  when  either  Det^  or  Det2  is  equal  to  zero.  As  in  the  cases 
for  isotropic  cylinders,  solutions  of  Detj  =  0  correspond  to  plane  strain  vibration,  and 
solutions  of  Det2  =  0  correspond  to  axial  shear  vibration.  When  the  variables  in  equa¬ 
tions  4.6  are  substituted  into  equation  4.9  using  normalized  elastic  constants,  we  find 
that  Det2  is  again  only  dependent  on  the  shear  wave  velocity  in  the  material. 

When  the  relationships  for  isotropic  materials,  equations  3.2,  are  substituted  into 
equations  4.6,  and  the  results  are  substituted  into  equation  4.9,  we  find  that  equation  4.9 
is  equivalent  to  equation  3.10.  Therefore,  the  normalized  cut-off  frequencies,  corre¬ 
sponding  to  axial  shear  modes,  for  a  transversely  isotropic  material  are  equal  to  the  nor¬ 
malized  cut-off  frequencies  for  an  isotropic  material.  This  means  that  the  normalized 
cut-off  frequencies  of  the  axial  shear  modes  of  vibration  in  the  cylinder  are  dependent 
only  on  the  shear  wave  velocity  in  the  material,  and  are  independent  of  any  fluid-load¬ 
ing  boundary  conditions,  at  least  when  the  cylinder  is  isotropic  or  transversely  isotro¬ 
pic. 

The  normalized  cut-off  frequencies  for  the  second,  third,  and  fourth  plane  strain 
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modes  of  the  fiber-reinforced  cylinder  were  computed  and  are  compared  to  the  normal¬ 
ized  cut-off  frequencies  of  the  soft  rubber  and  steel  cylinders  in  figure  11.  (Note  that 
figures  1 1-15  are  located  at  the  end  of  chapter  4,  starting  on  page  94).  The  cut-off  fre¬ 
quencies  are  normalized  to  the  shear  wave  velocity  in  each  material.  Using  the  material 
properties  for  rubber,  for  steel,  and  for  the  fiber-reinforced  material,  along  with  equa¬ 
tion  3.1,  the  shear  wave  velocities  in  the  three  materials  are 

Cj  =  72  m/sec,  for  rubber 
Cj  =  3218  m/sec,  for  steel 
Cj  =  75  m/sec,  for  fiber-reinforced  composite. 

In  general,  the  normalized  cut-off  fi^uencies  of  the  fiber-reinforced  cylinder  are 
very  similar  to  the  normalized  cut-off  frequencies  of  the  steel  cylinder.  For  the  second 
and  third  plane  strain  modes,  the  normalized  cut-off  frequencies  are  essentially  equal 
for  cylinders  with  a  ratio  of  inner  to  outer  radius  s  greater  than  or  equal  to  approximately 
0.7.  This  means  that,  at  zero  wavenumber  and  for  s  >  0.7 ,  the  hollow  fiber-reinforced 
cylinder  has  essentially  the  same  resonant  frequencies  for  these  plane  strain,  rigid  body 
modes  (i.e.  no  longitudinal  motion)  as  the  steel  cylinder. 

4.2.2  Frequency  spectrum 

The  frequency  spectrum  for  the  first  two  plane  strain  modes  and  the  first  axial  shear 
mode  of  the  fiber-reinforced  cylinder  with  a  ratio  of  the  inner  radius  to  the  outer  radius 
of  0.8,  and  for  real  wavenumbers  from  0  to  k,  is  shown  in  figure  12.  (Note  that  the  ap- 
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parent  piecewise  smooth  behavior  of  the  branches  in  figure  12  and  subsequent  figures 
is  due  to  conducting  the  root  searches  of  the  frequency  equation  at  discrete  wavenum¬ 
bers).  Again  we  have  labeled  the  branch  corresponding  to  the  first  plane  strain  mode  S 
the  branch  corresponding  to  the  first  axial  shear  mode  S2,  and  the  branch  corresponding 
to  the  second  plane  strain  mode  S3.  The  normalized  cut-off  frequency  for  the  S2  mode  is 


Q  =  0.223 


and  the  normalized  cut-off  frequency  for  the  S3  mode  is 

Q  =  0.5331. 

The  normalized  cut-off  frequency  for  the  S3  mode  is  about  1%  greater  than  the  normal¬ 
ized  cut-off  frequency  for  this  mode  in  the  steel  cylinder. 

At  low  wavenumbers,  the  normalized  phase  velocity  of  the  Sj  mode  is  approxi¬ 
mately  one-half  the  shear  wave  velocity  in  the  composite  material.  As  wavenumber  is 
increased,  the  phase  velocity  of  this  mode  asymptotically  approaches  the  shear  wave 
velocity  from  below. 

The  group  velocity  of  the  S3  mode  is  essentially  zero  below  5/n  =  0.01 .  Above 
this  wavenumber,  the  phase  velocity  of  this  mode  is  approximately  equal  to  the  longi¬ 
tudinal  wave  velocity  in  the  composite,  where  the  longitudinal  wave  velocity  is  defined 
as 
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(4.10) 


The  characteristics  of  the  S2  mode  are  quite  different  from  those  of  the  isotropic  ma¬ 
terials  studied  in  chapter  3.  At  low  wavenumbers,  this  mode  exhibits  weak  coupling  to 
the  S3  mode  and  the  normalized  phase  velocity  is  somewhat  greater  than  the  longitudi¬ 
nal  wave  velocity.  As  wavenumber  is  increased,  the  group  velocity  approaches  zero.  At 
wavenumbers  above  5/7C  =  0.1 ,  the  phase  velocity  of  this  mode  asymptotically  ap¬ 
proaches  the  shear  wave  velocity  from  above. 

Mirsky  (1965)  studied  the  axisymmetric  and  nonaxisymmetric  modes  of  wave 
propagation  in  transversely  isotropic  cylinders  and  compared  the  results  to  a  thick  shell 
theory.  In  his  analysis,  the  dynamic  behavior  of  beryllium  and  magnesium  was  com¬ 
pared  to  that  of  steel.  The  results  showed  that,  for  thin  cylinders,  ( s  =  0.9672 ),  the  fre¬ 
quency  spectra  for  the  steel,  beryllium  and  magnesium  cylinders  were  very  similar. 
However,  for  thick  cylinders,  (s  =  0.3333 ),  the  frequency  spectrum  for  the  beryllium 
cylinder  was  somewhat  different  from  that  of  the  steel  or  magnesium  cylinder.  It  should 
be  noted  that  in  the  low  frequency  region  of  the  frequency  spectrum  for  the  beryllium 
cylinder,  where  the  frequency  equation  contained  Bessel  functions  with  complex  argu¬ 
ments,  Mirsky  (1965)  substituted  numerical  results  from  a  thick  shell  theory  for  numer¬ 
ical  results  from  the  elasticity  theory. 


87 


To  the  author’s  knowledge,  a  comparison  of  the  dynamic  behavior  of  a  thick  rubber 
cylinder,  a  thick  steel  cylinder,  and  a  highly  anisotropic,  thick  fiber-reinforced  cylinder 
has  not  been  made  before  within  the  framework  of  the  theory  of  elasticity. 

4.3  A  HYPOTHETICAL  FLUID  COLUMN 

Next,  we  consider  the  cut-off  frequencies  and  dispersion  curves  of  a  hypothetical, 
cylindrical  fluid  column  with  pressure-release  (stress-free)  and  with  radial  displace¬ 
ment-free  boundary  conditions. 

4.3.1  Pressure-release  boundary  conditions 

The  frequency  equation  of  a  hypothetical  fluid  column,  within  the  transversely  iso¬ 
tropic  cylinder,  with  a  pressure-release  boundary  condition  is  identical  to  equation  3.13, 
that  is 


Jj  (aa)  =  0, 


(4.11) 


or. 


aa  =  3.832,7.016, 10.172, .... 


(4.12) 


4.3.2  Displacement-free  boundary  conditions 

The  frequency  equation  of  a  hypothetical  fluid  column,  within  the  transversely  iso- 
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tropic  cylinder,  with  a  displacement-free  boundary  condition  is  identical  to  equation 
3.16, that  is 


Jj  (aa)  -aaJ2(aa)  =  0, 


(4.13) 


or, 


aa  =  1.841,5.331,8.536,.... 


(4.14) 


Equations  4. 12  and  4. 14  represent  families  of  hyperbolas  with  finite  cut-off  fre¬ 
quencies.  The  first  two  cut-off  frequencies,  for  each  of  the  two  boundary  conditions  and 
s  =  0.8 ,  normalized  to  the  shear  wave  velocity  in  the  composite  material,  are 


=  19.08,  34.94,  pressure-release  boundary  conditions 
Q  =  9.17,  26.55,  displacement-free  boundary  conditions. 

These  normalized  cut-off  frequencies  are  very  similar  to  those  for  a  hypothetical  fluid 
column  in  the  rubber  cylinder. 

4.3.3  Frequency  spectrum 

The  first  two  branches,  corresponding  to  each  of  the  two  sets  of  boundary  condi- 
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tions  and  normalized  to  the  shear  wave  velocity  in  the  composite  material,  are  shown 
in  figure  13.  The  curves  labeled  1^  and  I2  represent  solutions  of  equations  4.12,  while 
the  curves  labeled  and  L2  represent  solutions  of  equations  4.14,  respectively. 

From  a  comparison  of  figures  12  and  13,  we  note  that  the  fluid  modes  exist  at  wave- 
numbers  and  frequencies  greater  than  the  solid  modes  of  the  cylinder.  This  means  that 
the  fluid  modes  wiU  not  couple  with  the  solid  modes,  at  least  in  the  wavenumber  and 
frequency  domains  of  interest.  This  is  similar  to  the  relationship  of  the  fluid  modes  and 
solid  modes  of  the  rubber  cylinder. 

4.4  A  FLUID-FILLED  CYLINDER 

The  next  case  to  be  investigated  is  that  of  a  fluid-filled  cylinder.  The  frequency 
equation  for  this  set  of  boundary  conditions  was  given  in  equation  3.18.  Again,  in  the 
limit  of  zero  wavenumber,  this  frequency  equation  separates  into  two  sub-determinants. 
The  solution  of  setting  one  sub-determinant  to  zero  corresponds  to  the  axial  shear 
modes,  which  as  stated  before,  are  found  to  be  dependent  only  on  the  shear  wave  ve¬ 
locity  in  the  material,  and  independent  of  any  fluid-loading  boundary  conditions. 

4.4.1  Cut-off  frequencies 

The  normalized  cut-off  frequencies  for  the  Sj,  S2,  and  S3  modes  for  the  fluid-filled 
fiber-reinforced  cylinder  are 


Q  =  0, 0.223,  and  0.398 , 


respectively. 

The  presence  of  the  internal  fluid  has  reduced  the  cut-off  frequency  of  the  S3  mode 
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by  approximately  25%.  No  additional  cut-off  frequencies,  corresponding  to  the  addi¬ 
tion  of  fluid  modes,  exist  within  the  normalized  frequency  band  of  0  to  7C.  The  degree 
of  reduction  of  the  cut-off  frequency  of  the  S3  mode,  and  the  absence  of  the  addition  of 
any  fluid  modes,  are  similar  to  the  characteristics  of  the  fluid-filled  rabber  cylinder. 

4.4.2  Frequency  spectrum 

The  frequency  spectrum  for  the  fluid-filled  fiber-reinforced  cylinder  is  shown  in 
figure  14.  From  a  comparison  of  figure  14  to  figure  12,  we  note  that  the  phase  velocity 
of  the  Sj  mode  of  the  fluid-filled  cylinder  is  less  than  the  phase  velocity  of  this  mode 
of  the  empty  cylinder,  for  all  wavenumbers.  For  example,  at  5  =  0.27t  the  phase  ve¬ 
locity  of  the  mode  is  reduced  approximately  35%  when  fluid  is  added  to  the  inside 
of  the  cylinder. 

The  phase  velocity  of  the  S3  mode  at  wavenumbers  greater  than  5/5T  =  0.01  is  vir¬ 
tually  unaffected  by  the  presence  of  the  fluid.  The  S2  mode  still  exhibits  weak  coupling 
to  the  S3  mode  at  low  wavenumbers.  The  phase  velocity  of  the  S2  mode  in  the  fluid- 
filled  cylinder  is  less  than  the  phase  velocity  of  this  mode  in  the  empty  cylinder  at  wave- 
numbers  between  0.02  <  5/7t  <  0.5 .  At  wavenumbers  greater  than  5/n  =  0.5  the 
phase  velocity  of  the  S2  mode  is  virtually  unchanged  and  asymptotically  approaches  the 
shear  wave  velocity  in  the  composite  material  from  above. 

4.5  A  CYLINDER  THAT  IS  FLUID-FILLED  AND  IMMERSED  IN  FLUID 

The  final  case  to  be  investigated  is  that  of  a  cylinder  that  is  fluid-filled  and  im¬ 
mersed  in  an  infinite  inviscid  fluid.  The  frequency  equation  for  this  set  of  boundary 
conditions  is  equation  3.22.  Again,  in  the  limit  of  zero  wavenumber,  this  frequency 
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equation  separates  into  two  sub-determinants  corresponding  to  the  axial  shear  modes 
and  plane  strain  modes.  The  solution  of  the  cut-off  frequencies  corresponding  to  the  ax¬ 
ial  shear  modes  remain  constant,  when  normalized  to  the  shear  wave  velocity  in  the  ma¬ 
terial,  and  independent  of  fluid-loading  boundary  conditions.  The  lowest,  non-zero  cut¬ 
off  frequency,  corresponding  to  a  plane  strain  mode,  is  the  cut-off  frequency  for  the  S3 
mode. 

4.5.1  Cut-off  frequencies 

The  normalized  cut-off  frequencies  for  the  Sj,  S2,  and  S3  modes  for  a  fiber-rein- 
forced  cylinder  that  is  fluid-filled  and  immersed  in  fluid  are 


Q  =  0, 0.223,  and  0.394. 


The  cut-off  frequency  of  the  S3  mode  is  essentially  unaffected  by  the  addition  of  the 
external  fluid  loading. 

4.5.2  Frequency  spectrum 

The  dispersion  curves  for  the  fiber-reinforced  cylinder  that  is  fluid-filled  and  im¬ 
mersed  in  fluid  are  shown  in  figure  15.  From  a  comparison  of  figures  14  and  15,  we 
note  that  the  S2  mode  is  essentially  unaffected  by  the  addition  of  the  external  fluid  load¬ 
ing.  The  phase  velocity  of  the  S^  mode  is  reduced  over  all  wavenumbers.  For  example, 
at  5/71  =  0.2  the  phase  velocity  has  been  reduced  by  approximately  18%. 

At  wavenumbers  less  than  5/n  =  0.05 ,  the  dispersion  characteristics  of  the  S3 
mode  are  unaffected  by  the  addition  of  the  external  fluid.  However,  at  higher  wavenum- 
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bers  the  phase  velocity  is  reduced  significantly.  For  example,  at  5/7i  =  0.1  the  phase 
velocity  is  reduced  by  46%  with  the  addition  of  the  external  fluid.  These  results  are  very 
different  from  the  results  for  the  rubber  cylinder  and  steel  cylinder. 

The  effect  of  the  external  fluid  on  the  dynamic  behavior  of  the  fluid-filled,  fiber- 
reinforced  cylinder  is  very  different  than  on  the  dynamic  behavior  of  the  fluid-filled 
rubber  and  steel  cylinder.  For  instance,  all  of  the  cut-off  frequencies  and  branches  of 
the  frequency  spectrum  for  the  fiber-reinforced  cylinder  are  real.  That  is,  unlike  the  re¬ 
sults  for  the  fluid-filled  steel  cylinder,  the  addition  of  the  external  fluid  has  not  intro¬ 
duced  any  damping  in  the  coupled  system.  Also,  unlike  the  results  for  the  fluid-filled 
rubber  cylinder,  the  addition  of  an  external  fluid  resulted  in  significant  reductions  in  the 
phase  velocities  for  the  Si  and  S3  modes. 

As  a  final  comment,  and  to  the  author’s  knowledge,  the  results  shown  in  figures  14 
and  15  represent  the  first  time  the  cut-off  frequencies  and  frequency  spectra  of  the  n=l, 
nonaxisymmetric  modes  of  wave  propagation  in  a  fluid-filled  cylinder  and  for  a  fluid- 
filled  cylinder  immersed  in  fluid,  have  been  computed  for  a  homogeneous,  transversely 
isotropic  cylinder. 
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Figure  11.  Comparison  of  the  cut-off  frequencies  of  the  plane  strain  modes  in  a  rubber, 
steel,  and  transversely  isotropic  cylinder  as  functions  of  the  ratio  of  inner  to  outer  radius 
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Figure  12.  Frequency  spectrum  of  the  transversely  cylinder  in  a  vacuum 
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Figure  13.  Frequency  spectrum  of  the  hypothetical  fluid  column  normalized  to  the  shear 
wave  velocity  in  the  composite  material 
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Figure  14.  Frequency  spectrum  of  the  fluid-filled  transversely  isotropic  cylinder 
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Figure  15.  Frequency  spectrum  of  the  transversely  isotropic  cylinder  that  is  fluid-filled 
and  immersed  in  fluid 
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5.  SUMMARY  AND  CONCLUSIONS 

The  frequency  equation  of  wave  propagation  in  a  transversely  isotropic  cylinder, 
that  is  fluid-filled  and  immersed  in  fluid,  was  developed  from  the  three-dimensional 
equations  of  elasticity  and  the  assumptions  of  perfect-slip  boundary  conditions  at  the 
fluid-solid  interfaces.  This  equation  is  general  in  axial  wavenumber  k,  circumferential 
wavenumber  n,  wall  thickness  h,  and  radian  frequency  co. 

The  frequency  equation  was  verified  by  specializing  the  variables  for  the  cases  of 
isotropic  cylinders  and  n  =  1  nonaxisymmetric  modes  and  comparing  the  normalized 
cut-off  frequencies  and  frequency  spectra  to  previous  results  for  hollow  cylinders  in  a 
vacuum  and  for  fluid-filled  cylinders. 

The  frequency  spectrum  for  the  n  =  1  modes  in  fluid-filled,  isotropic  cylinders  im¬ 
mersed  in  fluid  were  computed  and  compared  to  previous  results  for  similar  geometries 
and  for  the  n  =  0  modes.  It  was  found  that  the  effects  of  the  presence  of  the  inner  and 
outer  fluids  on  the  steel  cylinder,  for  the  n  =  1  modes,  were  similar  to  the  effects  for  the 
n  =  0  modes.  That  is,  the  inner  fluid  introduced  additional  modes  in  the  fluid-filled  cyl¬ 
inder  that  were  not  present  in  the  cylinder  in  a  vacuum.  The  outer  fluid  introduced 
damping  into  the  system  as  evidenced  by  the  presence  of  the  imaginary  term  in  the  com¬ 
plex  frequency. 

The  effects  of  the  inner  and  outer  fluids  on  the  soft  rubber  cylinder,  for  the  n  =  1 
modes,  were  quite  different  from  the  effects  on  the  steel  cylinder.  For  instance,  the  inner 
fluid  reduced  the  normalized  cut-off  frequencies  and  phase  velocities  of  the  plane  strain 
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modes  in  the  soft  rubber  cylinder  significantly  more  than  in  the  steel  cyhnder.  However, 
the  inner  fluid  did  not  introduce  any  additional  modes  within  the  wavenumber  and  fre¬ 
quency  domains  of  interest.  The  external  fluid  did  not  have  any  significant  additional 
effect  on  normalized  cut-off  frequencies  or  phase  velocities,  and,  unlike  the  case  for  the 
steel  cylinder,  did  not  introduce  any  damping  into  the  system. 

To  the  author’s  knowledge,  the  cut-off  frequencies  and  frequency  spectra  for  the 
n  =  1  modes  in  isotropic  cylinders,  with  fluid  loading  on  the  inner  and  on  outer  sur¬ 
faces,  have  not  been  computed  before  within  the  framework  of  the  theory  of  elasticity 

Next,  the  frequency  equation  was  specialized  for  the  cases  of  a  fiber-reinforced, 
transversely  isotropic  cylinder.  When  the  relationships  among  the  elastic  constants  of 
the  fiber-reinforced  cylinder  were  compared  to  those  of  some  naturally  occurring  hex¬ 
agonal  crystals,  the  fiber-reinforced  cylinder  was  found  to  be  highly  anisotropic. 

New  results  for  nonaxisymmetric  (n  =  1)  wave  propagation  in  a  transversely  isotro¬ 
pic  cylinder  were  obtained  for  the  cases  of  a  fluid-filled  cyhnder  and  a  cylinder  that  is 
fluid-filled  and  immersed  in  fluid.  The  normalized  cut-off  frequencies  and  frequency 
spectra  for  these  geometries  were  computed  and  compared  to  those  for  the  soft  rubber 
cylinder  and  the  steel  cylinder.  It  was  found  that  the  effects  of  the  inner  fluid  on  the 
wave  propagation  characteristics  of  the  fiber-reinforced  cylinder  were  similar  to  the  ef¬ 
fects  on  the  fluid-filled  rubber  cylinder,  as  evidenced  by  the  reductions  in  the  normal¬ 
ized  cut-off  frequencies  and  in  the  phase  velocities  of  the  plane  strain  modes.  However, 
the  effects  of  the  outer  fluid  on  the  two  types  of  cylinders  were  different.  Contrary  to 
the  case  of  the  soft  rubber  cyhnder,  the  presence  of  the  outer  fluid  on  the  fluid-filled. 
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fiber-reinforced  cylinder  further  reduced  the  phase  velocities  of  the  plane  strain  modes 
by  a  significant  amount. 

It  was  noted  that  normalized  cut-off  fi-equencies  corresponding  to  axial  shear  modes 
are  dependent  only  on  the  shear  wave  velocity  in  the  material  and  are  independent  of 
any  internal  or  external  fluid  loading  condition.  This  means,  at  zero  wavenumber  (cor¬ 
responding  to  rigid  body  motion),  the  normalized  resonant  frequency  for  harmonic  lon¬ 
gitudinal  vibration  is  the  same  for  isotropic  and  transversely  isotropic  cylinders, 
regardless  of  whether  they  are  fluid  filled  or  immersed  in  fluid  or  both. 

In  conclusion,  the  characteristics  of  wave  propagation  in  highly  anisotropic  cylin¬ 
ders,  with  or  without  fluid  loading,  may  be  markedly  different  from  the  characteristics 
of  wave  propagation  in  isotropic  cylinders.  Predictions  of  the  normalized  cut-off  fre¬ 
quencies  and,  in  particular,  phase  velocities  of  various  modes  in  a  transversely  isotropic 
cylinder,  based  on  results  for  isotropic  cylinders,  can  be  seriously  inaccurate. 
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6.  RECOMMENDATIONS  FOR  CONTINUED  INVESTIGATIONS 


1.  The  flexural  mode  is  often  the  mode  of  greatest  interest  in  nonaxisymmetric  wave 
propagation.  A  simplified,  low  wavenumber  model  for  this  mode  should  be  developed 
using  small  argument  approximations  for  the  Bessel  functions  in  the  frequency  equa¬ 
tion.  This  would  provide  a  method  to  quickly  analyze  the  effects  of  changes  in  material 
properties  of  the  transversely  isotropic  material  on  the  phase  velocity  of  this  mode. 

2.  The  numerical  computations  performed  in  this  dissertation  were  for  a  relatively 
thick  cylinder.  A  similar  study  of  the  nonaxisymmetric  wave  propagation  characteris¬ 
tics  in  fluid  loaded,  thin,  transversely  cylinders  should  also  be  conducted  and  compared 
to  thin  shell  theory. 

3.  Displacement  and  stress  distributions  in  the  fluids  and  the  cylinder  wall  should  be 
computed  over  a  range  of  wavenumbers  and  frequencies  in  order  to  analyze  the  degree 
of  coupling  of  the  radial,  circumferential  and  longitudinal  motions  in  the  cylinder  and 
the  pressure  distributions  in  the  inner  and  outer  fluids. 

4.  The  matrix  material  of  the  fiber-reinforced  cyhnder  was  assumed  to  be  linearly 
elastic.  Fiber-reinforced  materials  are  often  composed  of  a  linearly  viscoelastic  matrix 
material.  The  frequency  equation  for  the  fluid-loaded  cylinder  should  be  expanded  to 
include  the  effects  of  linear  viscoelasticity. 
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APPENDIX  A.  BOUNDARY  EQUATIONS 
The  boundary  equations,  which  were  formulated  in  chapter  2,  are 


[L]  {c}  =  0. 


(A.1) 


Equations  A.  1  constitute  a  system  of  eight,  linear,  homogeneous  algebraic  equations  in 
the  wave  amplitudes  {c}  ,  where 


{c}  —  ( Aj,  Bj,  A2,  B2,  Aj,  Bj,  Dj,  ©2)  , 


(A.2) 


and  where  the  components  of  [L]  are  given  below.  The  Yj's  are  +1  or  -1,  depending 
upon  the  type  of  Bessel  function  that  Wn  represents.  This  bookkeeping  variable  is  need¬ 
ed  due  to  the  differing  recursion  formulas  for  the  derivatives  of  the  Bessel  functions. 
For  instance,  Yj  =  -1 ,  when  represents  or  ,  and  Yi  =  1  when  Wn  rep¬ 

resents 


Lii  =  [2c66n(n-l) +[c44k^-pa)^-ic44Tlikja^W„(^ia) 

“Yi^Cgg^iaWjj^  j  (^1^) 

Li2  =  [2c66n(n-l) +[c44k^-p©^-ic44TiikJa^]z„(^ia) 

Li3  =  [2c66n(n-l) +[c44k^-p©^-ic44Ti2kJa^]w„(^2a) 
“'Y22C66^2^^n+  1  (^2^) 

Li4  =  [2C66n(n-l)  +[c44k^-p©^-ic44q2kja^]Zn(^2a) 
■*’2C66^2a^n+l  (^2^) 


(A.3) 


(A.4) 


(A.5) 


(A.6) 
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L,5  =  2c«n[{n-l)W„(43a)  +y3l3aW„^,  (Ija)] 

(A.7) 

1-16  =  2c,,n[(n-l)ZJ43a)-?3aZ„^,(43a)] 

(A.8) 

L,7  =  Pj<oVF„(^3a) 

(A.9) 

Li8  =  0 

(A.10) 

L21  =  (Tli  +  ik)a[nW„(^ja) +Yi^iaW„^j(^ja)] 

(A.11) 

L22  “  (^1  ■*■  ^k)  a  [nZjj  (^1^)  —  j  (^^a)  ] 

(A.12) 

L23  =  (Ti2  +  ik)a[nW^(^2^)  +  Y2^2^W^  +  i  (^2^)  1 

(A.13) 

L<24  ~  (112 ^  n  ^^2^^  ~  ^2^^n+ 1  ^^2^^  ^ 

(A.14) 

L25  =  inkaWj^(^3a) 

(A.15) 

L2g  =  inkaZ„(^3a) 

(A.  16) 

L27  =  0 

(A.17) 

^28  =  0 

(A.18) 

L3^  =  -2n[(n-l)WJ^ia)  +  j  (^^a)  ] 

(A.19) 

L32  =  -2n[(n-l)Z^(^^a)  -^iaZj^^j(^ja)] 

(A.20) 

L33  =  -2n[(n-l)W„(^2a)  +Y2^23W„  +  i(^2a)] 

(A.21) 

L34  =  -2n  [  (n-  1)  Z„  (^2^)  -  42aZ„.i  1 

(A.22) 

L35  =  [qV-2n(n-l)  jw„(^3a) 

(A.23) 

L36  =  (qV-2n(n-l)  jz„{^3a) -2^3aZ„^,  (^ja) 

(A.24) 

L37  =  0 

(A.25) 

^38  ~  ® 

(A.26) 

L4,  =nW„(^ia)-6Y,^iaW„^,(^,a) 

(A.27) 

L42  =  nZjj(^ja)  -^jaZjj_j_j  (^ja) 

(A.28) 
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^43  =  (^2^)  +  +  1  (^2^) 

L44  =  nZ„(^2a)-^2aZn.i(^2a)  (A-30) 

L45  =  i^W„(^3a)  (A.31) 

L46  =  *^ZJ^3a)  (A.32) 

L47  =  -(nW,(44a)  +  Y4^4aW„,i  (U))  (A.33) 

L43  =  0  (A.34) 

L51  =  [2c66n(n-l)  +[c44k^-p0)^-ic44^iky]w„(^ib) 

L52  =  [2c66n(n-l)  +(c44k^-p®^-iC44Tliky]Zn(^ib) 

''■^^66^1*’Zn+i  (^ib) 

L53  =  [2c66n(n-l)  +(c44k^-P®^-ic44Tl2kjb^W„(42t>)  ^^37) 

-Y22c,,^2bW„^i(^2b) 

L54  =  [2C66n(n-l)  +(c44k^-P®^-ic44Tl2ky]z„(^2b)  (^3g) 

■*■  + 1  (^2^^ 

L55  =  2Cjs''1("-1)W„(^3'>)  +Y353l>W„^l(l3b)]  (A.39) 

1-56  =  2‘:66'’[("-»Z„(^3‘>)-^3>>Zntl(^3'’» 

Lj,  =  0  (A.41) 

Ljs  =  P20>'‘b"G„(^5b)  (A.42) 

L„  =  (nl-^ik)b[nW„(5,b)-^Y,^,bW„^l(^,b)]  (A.43) 

L6j=  (Tl,+ik)b[nZ„(4ib)-4,bZ„3.,(4,b)]  (A.44) 

=  (Hj  -Hk)  b  [nW  „  (Ijb)  +  ,  (4jb)  ]  (A.45) 
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L„  =  (Tlj  +  ilc)b[nZ„(|2b)-5jbZ„^j(^jb)  ] 

(A.46) 

Lj5  =  inlcbW„(^3b) 

(A.47) 

^66  “  (^2^) 

(A.48) 

^67  ~  ^ 

(A.49) 

^68  “  ^ 

(A.50) 

L71  =  -2n[(n-l)W„(^ib)  +Yi^ibW„^i(^ib)] 

(A.51) 

L72  =  -2n[(n-l)Z^(^^b)  -^jbZ^^j(^,b)] 

(A.52) 

L73  =  -2n[(n-l)W^(^2b)  +Y2^2bW„^i(^2b)] 

(A.53) 

L74  =  -2n[(n-l)Z^(^2b)  -^2bZn^.i  (^2^)1 

(A.54) 

L,5  =  (qV-2n(n-l)  Jw^C^jb)  +2Y343bW„^,  (Ijb) 

(A.55) 

L76  =  [qV-2n(n-l)  Jz„(53b) -2|3bZ„^,  (^jb) 

(A.56) 

L77  =  0 

(A.57) 

^78  ~  ® 

(A.58) 

L«l  ="W„(|,b)+Y,^,bW„^,(?,b) 

(A.59) 

L83  =  nZ„(^jb)-^,bZ„,,(^,b) 

(A.60) 

Lgs  =  nWJIjb)  +Y2l2bW.„(42b) 

(A.61) 

L34  =  nZ„(^3b)-42bZ„^,(?3b) 

(A.62) 

Lss  =  ”W„  (^jb) 

(A.63) 

L86  =  "W„  (43b) 

(A.64) 

Lg?  =  0 

(A.65) 

Lss  =  -(-WJ^sb)  +Y5^5bW„^,(45b)) 

(A.66) 
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The  solutions  of  A.  1  are  nontrivial  when  the  determinant  of  the  coefficients  of  the 
wave  amplitudes  {c}  vanishes,  i.e.,  when 


|L|  =  0. 


(A.67) 


Equation  A.67  is  a  complex  algebraic  equation.  The  real  and  imaginary  parts  must  van¬ 
ish  simultaneously.  To  aid  in  the  computation  of  the  solutions  of  equation  A.67 ,  all  the 
physical  variables  are  normalized  as  follows. 


a 

b 

h 


h 


s 

(1-s) 


h— - T  ,  where 

(1-s) 

(b-a),and 


s  = 


a 

b’ 


(A.68) 


and 


SI 

5 

c 


=  — ,  normalized  fi’equency, 

Ct 

=  kh,  normalized  wavenumber, 

=  ^  ,  normalized  phase  speed, 
o 


(A.69) 


The  branches  of  the  frequency  spectrum  are  calculated  by  selecting  a  real  value  of 
the  normalized  wavenumber  and  searching  for  the  real  or  complex  value  of  the  normal¬ 
ized  frequency.  A  Mathematica  (Wolfram,  1991)  program  was  developed  to  conduct 
the  root  searches  using  Mueller’s  method  (Press  et  al,  1986).  Mueller’s  method  is  an 
algorithm  that  finds  a  real  or  complex  zero  of  a  complex  function,  which  in  this  case  is 
equation  A.67.  An  example  Mathematica  program  is  provided  in  appendix  B. 
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APPENDIX  B.  MATHEMATICA  PROGRAM  TO  CONDUCT  ROOT 


SEARCHES  OF  THE  FREQUENCY  EQUATION 


This  is  a  Mathematica  script  file  to  conduct  root  searches  using  Mueller’s  method.  The 
comments  and  instructions  are  in  Italic.  The  Mathematica  commands  can  be  typed  or 
cut  and  pasted  into  a  Mathematica  session. 

The  following  define  the  ratio  of  inner  to  outer  radius,  the  circumferential  wavenumber, 
the  elastic  constants  and  density  of  the  cylinder,  and  the  densities  and  acoustic  phase 
velocities  of  the  inner  and  outer  fluids. 

s=.8;n=l; 
cl  1  =3.7  c44; 
cl2  =  1.7  c44; 
cl3  =  2.3  c44; 
c33  =  100  c44; 
c66  =  c44; 
rho  =  1014; 
rhol =  1000; 
rho2  =  1000; 
cl = 1500; 
c2  =  1500; 

The  following  Mathematica  instructions  simplify  the  arguments  of  the  Bessel  functions 
and  their  coefficients. 

kappal  =  l/2((cll+c44)/cll  omega^2  +  (cl3''2  +  2  cl3  c44  -  cll  c33)/(cll  c44)  * 
delta''2); 

kappal  =  Cancel[kappal]; 
kappal  =  Simplify[kappal]; 

kappa2  =  (c44/cll  omega''2  -  c33/cll  delta^2)(omega^2  -  delta^2); 
kappa3  =  Simplify  [kappal  -  Sqrt[kappal''2  -  kappa2]]; 
kappa3  =  PowerExpand[kappa3]; 
kappa3  =  Simplify[kappa3]; 

kappa4  =  Simplify  [kappal  +  Sqrt[kappal''2  -  kappa2]]; 

kappa4  =  PowerExpand[kappa4]; 

kappa4  =  Simplify[kappa4]; 

kappa5  =  c44/c66  (omega'^2  -  delta''2); 
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kappa6  =  -(cl3+c44)  delta''2  kappa3/(c44  omega'^2  -  c33  delta'^2  -  c44  kappa3); 
kappa?  =  (c44  delta''2  -  c44  omega'^2  +  cll  kappa4)/(cl3  +  c44); 


denom  =  c44  omega''2-c33  delta^2  -c44  kappa3; 
denom  =  Collect[denom,c44]; 
denom  =  Simplify  [denom]; 

numl  =  (c33+c44)  delta^2  omega''2  -  c33  deltaM  -  c44  omegaM  + 

(c44  omega^2  +  cl3  delta'^2)  kappa3; 
numl  =  Collect[numl,c44]; 
numl  =  Simplify[numl]; 
numl=  Cancel[numl/denom]; 

nuni2  =  (cl3(delta''2  -  omega''2)  -cll  kappa4)/(cl3  +  c44); 
num2  =  Cancel  [num2]; 
num2  =  Simplify  [num2]; 

num3  =  delta''2(c33  delta''2  -  c44  omega'^2  -  cl3  kappa3); 
num3  =  Collect[num3,c44]; 
num3  =  Cancel  [num3/denom]; 

num4  =  (cl  1  kappa4  -  c44  omega''2  -  cl3  delta^2)/(cl3+c44); 
num4  =  Cancel[num4]; 
num4  =  Simplify[num4]; 
c44  =  5.7*10^6; 

alpha  =  c44/(rho*cl''2)  omega^2  -  delta''2; 
beta  =  c44/(rho*c2'^2)  omega'^2  -  delta'^2; 

xil  =  Sqrt[kappa3]; 
xi2  =  Sqrt[kappa4]; 
xi3  =  Sqrt[kappa5]; 
xi4  =  Sqrt[Abs[alpha]]; 
xi5  =  Sqrt[Abs[beta]]; 

The  next  Mathematica  instructions  define  the  components  of  the  determinant  ILj. 

Lll  =  (2c66/c44  n  (n-1)  +  numl  (s/(l-s))^2)* 

BesseU[n,xil  s/(l-s)]  +  2c66/c44  xil  s/(l-s)  BesseU[n+l,xil  s/(l-s)]; 

L12  =  (2c66/c44  n  (n-1)  +  numl  (s/(l-s))^2)* 

BesselY[n,xil  s/(l-s)]  +  2c66/c44  xil  s/(l-s)  Bessel Y[n+1, xil  s/(l-s)]; 

L13  =  ((2c66/c44  n  (n-1)  +  num2  (s/(l-s))'^2)* 

BesseU[n,xi2  s/(l-s)]  +  2c66/c44  xi2  s/(l-s)  BesselJ[n+l,xi2  s/(l-s)]); 

L14  =  ((2c66/c44  n  (n-1)  -h  num2  (s/(l-s))''2)* 
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BesselY[n,xi2  s/(l-s)]  +  2c66/c44  xi2  s/(l-s)  BesselY[n+l,xi2  s/(l-s)]); 

L15  =  2c66/c44  n(  (n-1)  BesselJ[n,xi3  s/(l-s)]  -  xi3  s/(l-s)* 
BesseU[n+l,xi3  s/(l-s)]); 

L16  =  2c66/c44  n(  (n-1)  BesselY[n,xi3  s/(l-s)]  -  xi3  s/(l-s)* 
BesselY[n+l,xi3*s/(l-s)]); 

L17  =  (rhol/rho  omega'^2)  (s/(l-s))'^2  Bessell[n,xi4  s/(l-s)]; 


L18  =  0; 


L21  =  num3  (n  BesselJ[n,xil  s/(l-s)]  -  xil  s/(l-s)  BesselJ[ 
n-t-1,  xil  s/(l-s)]); 

L22  =  num3  (n  BesselY[n,xil  s/(l-s)]  -  xil  s/(l-s)  BesselY[ 
n+1,  xil  s/(l-s)]); 

L23  =  num4  (n  Bessel J[n,xi2  s/(l-s)]  -  xi2  s/(l-s)  BesseU[ 
n-(-l,  xi2  s/(l-s)]); 

L24  =  num4  (n  BesselY[n,xi2  s/(l-s)]  -  xi2  s/(l-s)  BesselY[ 
n+1,  xi2  s/(l-s)]); 

L25  =  -n  delta^2  BesselJ[n,xi3  s/(l-s)]; 

L26  =  -n  delta''2  BesselY[n,xi3  s/(l-s)]; 


L27  =  0; 


L28  =  0; 


L31  =  c66/c44  (-2n((n-l)  BesseU[n,xil  s/(l-s)]  - 
xil  s/(l-s)  Bessel J[n+1, xil  s/(l-s)])); 

L32  =  c66/c44  (-2n((n-l)  BesselY[n,xil  s/(l-s)]  - 
xil  s/(l-s)  BesselY[n+l,xil  s/(l-s)])); 

L33  =  c66/c44  (-2n((n-l)  BesseU[n,xi2  s/(l-s)]  - 
xi2  s/(l-s)  BesselJ[n+l,xi2  s/(l-s)])); 

L34  =  c66/c44  (-2n((n-l)  BesselY[n,xi2  s/(l-s)]  - 
xi2  s/(l-s)  Bessel Y[n+l,xi2  s/(l-s)])); 
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L35  =  c66/c44  ((kappa5  (s/(l-s))'^2  -  2n  (n-1))  BesselJ[n,xi3* 
s/(l-s)]  -2  xi3  s/O-s)  BesselJ[n+l,xi3  s/(l-s)]); 

L36  =  c66/c44  ((kappaS  (s/(l-s))''2  -  2n  (n-1))  BesselY[n,xi3* 
s/(l-s)]  -2  xi3  ^h-s)  BesselY[n+l,xi3  s/(l-s)]); 


L37  =  0; 


L38  =  0; 


L41  =  n  BesselJ[n,xil  s/(l-s)]  -  xil  s/(l-s)  Bessel J[n-i-l,xil  s/(l-s)]; 

L42  =  n  BesselY[n,xil  s/(l-s)]  -  xil  s/(l-s)  Bessel Y[n+1, xil  s/(l-s)]; 
L43  =  n  BesselJ[n,xi2  s/(l-s)]  -  xi2  s/(l-s)  BesseU[n+l,xi2  s/(l-s)]; 

L44  =  n  BesselY[n,xi2  s/(l-s)]  -  xi2  s/(l-s)  BesselY[n+l,xi2  s/(l-s)]; 
L45  =  n  BesselJ[n,xi3  s/(l-s)]; 

L46  =  n  Bessel Y[n,  xi3  s/(l-s)]; 

L47  =  -(n  Bessell[n,xi4  s/(l-s)]  +  xi4  s/(l-s)  Bessell[n+l,xi4  s/(l-s)]); 
L48  =  0; 

L51  =  (2c66/c44  n  (n-1)  +  numl  (l/(l-s))^2)* 

BesseU[n,xil  l/(l-s)]  +  2c66/c44  xil  l/(l-s)  Bessel J[n+1, xil  l/(l-s)]; 

L52  =  (2c66/c44  n  (n-1)  +  numl  (l/(l-s))'^2)* 

BesselY[n,xil  l/(l-s)]  +  2c66/c44  xil  l/(l-s)  Bessel Y[n+1, xil  l/(l-s)]; 

L53  =  ((2c66/c44  n  (n-1)  +  nuin2  (l/(l-s))^2)* 

BesselJ[n,xi2  l/(l-s)]  +  2c66/c44  xi2  l/(l-s)  BesseU[n+l,xi2  l/(l-s)]); 

L54  =  ((2c66/c44  n  (n-1)  +  num2  (l/(l-s))^2)* 

BesselY[n,xi2  l/(l-s)]  +2c66/c44  xi2  l/(l-s)  Bessel Y[n+l,xi2  l/(l-s)]); 

L55  =  2c66/c44  n(  (n-1)  BesselJ[n,xi3  l/(l-s)]  -  xi3  l/(l-s)  * 
BesseU[n+l,xi3  l/(l-s)]); 

L56  =  2c66/c44  n(  (n-1)  BesselY[n,xi3  l/(l-s)]  -  xi3  l/(l-s)  * 
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BesselY[n+l,xi3  l/(l-s)]); 

L57  =  0; 

L58  =  (rho2/rho)  omega^2  (l/(l-s))'^2  BesselK[n,xi5  l/(l-s)]; 

L61  =  num3  (n  BesselJ[n,xil  l/(l-s)]  -  xil  l/(l-s)  BesselJ[ 
n+1,  xil  l/(l-s)]); 

L62  =  num3  (n  BesselY[n,xil  l/(l-s)]  -  xil  l/(l-s)  BesselY[ 
n+1,  xil  l/(l-s)]); 

L63  =  num4  (n  BesseU[n,xi2  l/(l-s)]  -  xi2  l/(l-s)  BesselJ[ 
n+1,  xi2  l/(l-s)]); 

L64  =  num4  (n  BesselY[n,xi2  l/(l-s)]  -  xi2  l/(l-s)  BesselY[ 
n+1,  xi2  l/(l-s)]); 

L65  =  -n  delta''2  BesselJ[n,xi3  l/(l-s)]; 

L66  =  -n  delta'^2  BesselY[n,xi3  l/(l-s)]; 

L67  =  0; 

L68  =  0; 

L71  =  c66/c44  (-2n((n-l)  BesselJ[n,xil  l/(l-s)]  -  xil  l/(l-s)* 
BesseU[n+l,xil  l/(l-s)])); 

L72  =  c66/c44  (-2n((n-l)  BesselY[n,xil  l/(l-s)]  -  xil  l/(l-s)  * 
BesselY[n+l,xil  l/(l-s)])); 

L73  =  c66/c44  (-2n((n-l)  BesselJ[n,xi2  l/(l-s)]  -  xi2  l/(l-s)* 
BesselJ[n+l,xi2  l/(l-s)])); 

L74  =  c66/c44  (-2n((n-l)  BesselY[n,xi2  l/(l-s)]  -  xi2  l/(l-s)* 
Bessel Y[n+l,xi2  l/(l-s)])); 

L75  =  c66/c44  ((kappa5  (l/(l-s))''2  -  2n  (n-1))  BesselJ[n,xi3* 
l/(l-s)]  -2  xi3  l/(l-s)  BesselJ[n+l,xi3  l/(l-s)]); 

L76  =  c66/c44  ((kappa5  (l/(l-s)y'2  -  2n  (n-1))  BesselY[n,xi3* 

l/(l-s)]  -2xi3  l/(l-s)  BesselY[n+l,xi3  l/(l-s)]); 
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L77  =  0; 


L78  =  0; 

L81  =  n  BesselJ[n,xil  l/(l-s)]  -  xil  l/(l-s)  BesselJ[n+l,xil  l/(l-s)]; 

L82  =  n  BesselY[n,xil  l/(l-s)]-xil  l/(l-s)  BesselY[n+l,xil  l/(l-s)]; 

L83  =  n  BesselJ[n,xi2  l/(l-s)]  -  xi2  l/(l-s)  BesseU[n+l,xi2  l/(l-s)]; 

L84  =  n  BesselY[n,xi2  l/(l-s)]  -  xi2  l/(l-s)  BesselY[n+l,xi2  l/(l-s)]; 

L85  =  n  BesselJ[n,xi3  l/(l-s)]; 

L86  =  n  BesselY[n,xi3  l/(l-s)]; 

L87  =  0; 

L88  =  -(n  BesselK[n,xi5  l/(l-s)]  -  xi5  l/(l-s)  BesselK[n+l,xi5  l/(l-s)]); 

My  root”  is  a  Mathematica  program  (somewhat  like  a  subroutine)  that  is  called  with 
the  arguments  x,  y,  z.  These  arguments  are  the  "first  guess”  root  locations.  This  pro¬ 
gram  uses  Mueller’s  method  to  locate  the  closest  root  to  the  first  guess.  If  a  root  is  not 
found  after  25  iterations,  the  program  stops.  The  criteria  for  convergence  is  that  the 
difference  between  the  current  value  of  a  possible  root  and  the  last  calculated  value  is 
less  than  1x10^10. 

myroot  :=  For[i=l,i<25,i++,( 

q  =  (z-y)/(y-x); 

omega  =  x; 

xil  =  Sqrt[kappa3]; 

xi2  =  Sqrt[kappa4]; 

xi3  =  Sqrt[kappa5]; 

xi4  =  Sqrt[Abs[alpha]]; 

xi5  =  Sqrt[Abs[beta]]; 

matrix  =  {{N[L11],N[L12],N[L13],N[L14],N[L15],N[L16],N[L17],N[L18]}, 
{N[L21],N[L22],N[L23],N[L24],N[L25],N[L26],N[L27],N[L28]}, 
{N[L31],N[L32],N[L33],N[L34],N[L35],N[L36],N[L37],N[L38] } , 
{N[L41],N[L42],N[L43],N[L44],N[L45],N[L46],N[L47],N[L48]}, 
{N[L51],N[L52],N[L53],N[L54],N[L55],N[L56],N[L57],N[L58] } , 
{N[L61],N[L62],N[L63],N[L64],N[L65],N[L66].N[L67],N[L68]}, 

{N[L7 1],N[L72],N[L73],N[L74],N[L75],N[L76],N[L77],N[L78] } , 
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{N[L81],N[L82],N[L83],N[L84],N[L85],N[L86],N[L87],N[L88] } }; 
px  =  Det[matrix]; 
omega  =  y; 
xil  =  Sqrt[kappa3]; 
xi2  =  Sqrt[kappa4]; 
xi3  =  Sqrt[kappa5]; 
xi4  =  Sqrt[Abs[alpha]]; 
xi5  =  Sqrt[Abs[beta]]; 

matrix  ={  {N[L11],N[L12],N[L13],N[L14],N[L15],N[L16],N[L17],N[L18]}, 
{ N[L21]  ,N[L22],N[L23]  ,N[L24],N[L25]  ,N[L26],N[L27]  ,N[L28] } , 
{N[L31],N[L32],N[L33],N[L34],N[L35],N[L36],N[L37],N[L38]}, 

{ N[L4 1]  ,N[L42],N[L43]  ,N[L44],N[L45]  ,N[L46],N[L47]  ,N[L48] } , 
{N[L51],N[L52],N[L53],N[L54],N[L55],N[L56],N[L57],N[L58]}, 
{N[L61],N[L62],N[L63],N[L64],N[L65],N[L66],N[L67],N[L68]}, 

{ N[L7 1]  ,N[L72],N[L73]  ,N[L74],N[L75]  ,N[L76],N[L77]  ,N[L78] } , 
{N[L81],N[L82],N[L83],N[L84],N[L85],N[L86],N[L87],N[L88] } }; 
py  =  Det[matrix]; 
omega  =  z; 
xil  =  Sqrt[kappa3]; 
xi2  =  Sqrt[kappa4]; 
xi3  =  Sqrt[kappa5]; 
xi4  =  Sqrt[Abs[alpha]]; 
xi5  =  Sqrt[Abs[beta]]; 

matrix  =  {{N[L11],N[L12],N[L13],N[L14],N[L15],N[L16],N[L17],N[L18]}, 
{ N[L2 1  ]  ,N[L22]  ,N[L23]  ,N[L24],N[L25]  ,N[L26]  ,N[L27]  ,N[L28] } , 
{N[L31],N[L32],N[L33],N[L34],N[L35],N[L36],N[L37],N[L38]}, 
{N[L41],N[L42],N[L43],N[L44],N[L45],N[L46],N[L47],N[L48]}, 
{N[L51],N[L52],N[L53],N[L54],N[L55],N[L56],N[L57],N[L58]}, 
{N[L61],N[L62],N[L63],N[L64],N[L65],N[L66],N[L67],N[L68] } , 
{N[L71],N[L72],N[L73],N[L74],N[L75],N[L76],N[L77],N[L78]}, 
{N[L81],N[L82],N[L83],N[L84],N[L85],N[L86],N[L87],N[L88] } }; 
pz  =  Det[matrix]; 
a  =  q  pz  -  q  (1+q)  py  +  q''2  px; 
b  =  (2q  +  1)  pz  -  (1  +  q)'^2  py  +  q''2  px; 
c  =  (l  +q)  pz; 
deni  =  b  +  Sqrt[b''2  -  4  a  c]; 
den2  =  b  -  Sqrt[b''2  -  4  a  c]; 
d  =  Abs[denl]; 
e  =  Abs[den2]; 

If[d  >  e,  denom  =  deni,  denom  =  den2]; 
next  =  Chop[z  -  (z  -  y)  (2  c/  denom),  10^-6]; 

X  =  y;  y  =  z;  z  =  next; 
f  =  Chop[z-y]; 
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If[f  !=  0,  ContinueO,  Retum[next]]; 

omega  =  next; 

xil  =  Sqrt[kappa3]; 

xi2  =  Sqrt[kappa4]; 

xi3  =  Sqrt[kappa5]; 

xi4  =  Sqrt[Abs[alpha]]; 

xi5  =  Sqrt[Abs[beta]]; 

matrix  =  {{N[L11],N[L12],N[L13],N[L14],N[L15],N[L16],N[L17],N[L18]}, 
{N[L21],N[L22],N[L23],N[L24],N[L25],N[L26],N[L27],N[L28]}, 
{N[L31],N[L32],N[L33],N[L34],N[L35],N[L36],N[L37],N[L38]}, 
{N[L41],N[L42],N[L43],N[L44],N[L45],N[L46],N[L47],N[L48] } , 
{N[L51],N[L52],N[L53],N[L54],N[L55],N[L56],N[L57],N[L58]}, 
{N[L61],N[L62],N[L63],N[L64],N[L65],N[L66],N[L67],N[L68]}, 
{N[L71],N[L72],N[L73],N[L74],N[L75],N[L76],N[L77],N[L78]}, 
{N[L81],N[L82],N[L83],N[L84],N[L85],N[L86],N[L87],N[L88] } }; 
g  =  Chop[Det[matrix]]; 

If[g  !=  0,Contmuen,Retum[next]])] 

The  following  is  a  representative  set  of  Mathematica  instructions  to  locate  two  roots  for 
two  branches.  Typically,  a  root  if  found  in  6  or  7  iterations.  Several  of  these  sets  of  in¬ 
structions  can  be  input  at  once  to  run  the  root  searches  unattended.  The  answers,  which 
are  displayed  on  the  screen,  and  can  be  saved  in  a  text  file,  are  shown  in  Italics. 

delta  =  N[.005*Pi] 
x=.26;y=.27;z=.28;myroot 
omega  =  0.2705293700078156 
x=.39;y=.395;z=.4;myroot 
omega  =  0.3950075764334804 

The  numerical  computations  were  performed  using  machine  precision  numbers,  where 
precision  r^ers  to  the  total  number  of  significant  decimal  digits.  In  the  case  of  the  com¬ 
puters  used  to  compute  the  results  shown  in  the  body  of  the  dissertation,  the  machine 
precision  was  16.  The  accuracy  of  a  root  is  the  number  of  significant  decimal  digits  to 
the  right  of  the  decimal  point. 
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APPENDIX  C.  MATERIAL  PROPERTIES 


Steel 


Rubber 


10  2 

Young's  modulus:  21  x  10  N/m 
Poisson's  ratio:  0.3 

Density:  7800  kg/m^ 


7  2 

Young's  modulus:  1.5  x  10  N/m 
Poisson's  ratio:  0.45 

Density:  1000  kg/m^ 


Inner  and  Outer  Fluids 


Acoustic  phase  velocity:  1500  m/sec 
Density:  1000  kg/m^ 


2 


Composite  material  matrix 


7  2 

Young's  modulus:  1.5  x  10  N/m 
Poisson's  ratio:  0.45 

Density:  1000  kg/m^ 
Volume  fraction:  0.9 
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I 


Composite  material  reinforcing  fibers 


9  2 

Young's  modulus:  5.5  x  10  N/m 
Poisson's  ratio:  0.15 

3 

Density:  1 140  kg/m 
Volume  fraction:  0.1 
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APPENDIX  D.  MATHEMATICA  PROGRAM  TO  DETERMINE 


THE  EFFECTIVE  ELASTIC  CONSTANTS 


This  is  a  Mathematica  script  file  to  determine  the  elastic  constants  of  a  homogeneous, 
transversely  isotropic  material.  The  comments  and  instructions  in  Italic  are  not  Math¬ 
ematica  commands.  The  Mathematica  commands  can  be  typed  or  cut  and  pasted  into  a 
Mathematica  session  to  produce  the  results  shown  in  Italics. 

First,  we  define  Kronecker’s  delta,  the  unit  vector  {m},  and  the  symmetry  of  the  com¬ 
ponents  of  the  strain  tensor. 

kroneckerDelta[i_,j_]  :=  lf[i==j,l,0] 

m[p_,qj  :=  m[p]  m[q] 

eps[2,l]=eps[l,2] 

eps[l,3]=eps[3,l] 

eps[3,2]=eps[2,3] 

The  following  is  equation  2.6. 

stress[i_j_]  :=  Expand[  kroneckerDelta[i,j]  (cl  Sum[eps[kjc],{k,l,3)]  + 

c2  Sum[m[p,q]  eps[p,q],{p,l,3},{q,l,3}])  +  ni[ij]  (c2  Sum[eps[k,k],{k,l,3)]  + 
c3  Sum[m[p,q]  eps[p,q],  {p,l,3},{q,l,3}])  +  c4  eps[io]  + 
c5  (Sum[m[i,p]  eps[p,j],  {p,l,3}]  +  Sum[m[j,p]  eps[p,i],{p,l,3)])  1; 

When  the  direction  cosines  are  defined  as  mil]  =  0,  m[2]  =  0,  m[3]  =  I,  the  elastic 
constants  are  found  from  the  equations  below. 

m[l]=0 
m[2]=0 
m[3]  =  1 


cll  =  Simplify[Coefficient[stress[l,l],eps[l,l]]] 
cl2  =  Simplify[Coefficient[stress[l,l],eps[2,2]]] 
cl3  =  Simplify[Coefficient[stress[l,l],eps[3,3]]] 
cl4  =  Simplify[Coefficient[stress[l,l],eps[2,3]]] 
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cl5  =  Simplify[Coefficient[stress[l,l],eps[3,l]]] 
cl6  =  Simplify[Coefficient[stress[l,l],eps[l,2]]] 

c21  =  Simplify[Coefficient[stress[2,2],eps[l,l]]] 
c22  =  Simplify[Coefficient[stress[2,2],eps[2,2]]] 
c23  =  Simplify[Coefficient[stress[2,2],eps[3,3]]] 
c24  =  Simplify[Coefficient[stress[2,2],eps[2,3]]] 
c25  =  Simplify[Coefficient[stress[2,2],eps[3,l]]] 
c26  =  Simplify[Coefficient[stress[2,2],eps[l,2]]] 

c31  =  Simplify[Coefficient[stress[3,3],eps[l,l]]] 
c32  =  Simplify  [Coefficient[stress[3,3].eps[2,2]]] 

c33  =  Simplify[Coefficient[stress[3,3],eps[3,3]]] 
c34  =  Simplify[Coefficient[stress[3,3].eps[2,3]]] 
c35  =  Simplify[Coefficient[stress[3,3],eps[3,l]]] 

c36  =  Simplify[Coefficient[stress[3,3]»eps[l,2]]] 

c41  =  Simplify  [Coefficient[stress[2, 3]  ,eps[  1,1]]] 
c42  =  Simplify[Coefficient[stress[2,3],eps[2,2]]] 
c43  =  Simplify[Coefficient[stress[2,3],eps[3,3]]] 
c44  =  Simplify  [Coefficient[stress[2, 3]  ,eps[2,3]]/2] 
c45  =  Simplify[Coefficient[stress[2,3],eps[3,l]]/2] 
c46  =  Simplify[Coefficient[stress[2,3],eps[l,2]]/2] 

c51  =  Simplify[Coefficient[stress[3,l],eps[l,l]]] 
c52  =  Simplify[Coefficient[stress[3,l],eps[2,2]]] 
c53  =  Simplify[Coefficient[stress[3,l],eps[3,3]]] 
c54  =  Simplify[Coefficient[stress[3,l],eps[2,3]]/2] 
c55  =  Simplify[Coefficient[stress[3,l],eps[3,l]]/2] 
c56  =  Simplify[Coefficient[stress[3,l],eps[l,2]]/2] 

c61  =  Simplify[Coefficient[stress[l,2],eps[l,l]]] 
c62  =  Simplify[Coefficient[stress[l,2],eps[2,2]]] 
c63  =  Simplify  [Coefficient[stress[l, 2]  ,eps[3,3]]] 
c64  =  Simplify[Coefficient[stress[l,2],eps[2,3]]/2] 
c65  =  Simplify[Coefficient[stress[l,2],eps[3,l]]/2] 
c66  =  Simplify[Coefficient[stress[l,2],eps[l,2]]/2] 

The  answers  are  as  follows: 

cll  =  cl  +  c4 
cl2  =  cl 
cl3  =  cl  +  c2 


119 


c33  =  cl  +  2c2 
c44  =  (c4  +  c5)l2 
c66  =  c4l2 

When  the  following  relationships  between  the  constants  cl,  c2,  c3,  c4,  and  c5  and  the 
more  familiar  engineering  constants  are  defined, 

cl  =bt  -  gt 

c2  =  (2  vl  -  l)bt  +  gt 

c3  =E1  +  (2  vl  - 1)'^2  bt-  4  gl  +  gt 

c4  =  2  gt 

c5  =  2  gl  -  2  gt 

the  elastic  constants,  in  terms  of  the  engineering  constants  simplify  to: 

cll  =  bt  +  gt 
cl 2  =  bt-  gt 
cl3  =  2btvl 
c33  =  El  +  4  bt  vl^2 
c44  =  gl 
c66  =  gt 

When  the  cylinder  is  made  from  a  fiber  reinforced  matrix, 
the  ^ective  properties  can  be  found  from 

El  =  Ef  Vf  +  Em  Vm 
Et  =  (Ef  Em)/(Vm  Ef  +  Vf  Em) 
vl  =  Vm  vm  +  Vf  vf 
vt  =  vm 

gl  =  (gm  gf)/(Vm  gf  +  Vf  gm) 
gt  =  Et/(2(l+vt)) 
rho  =  pf  Vf  +  pm  Vm 
bt  =  .5  El/(l-vt-2vl'^2  Et/El) 

The  following  properties  were  used  for  a  soft  rubber  matrix 
reinforced  with  nylon  cords 

Ef  =  5.5*10^9 
Em  =  1.5*10^7 
vm  =  0.45 
vf  =  0.15 

gm  =  Em/(2(l+vm)) 
gf  =  Ef/(2(l+vf)) 
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pm  =  1000 
pf=  1140 
Vf  =  0.10 
Vm  =  0.90 

The  resultant  effective  properties  are 

El  =  5.635*10^ 

Et  =  1. 66616*  m? 
vl  =  0.42 
vt  =  0.45 

gl  =  5.74575*10'^ 
gt  =  5.74539*10^ 
bt  =  154398*10^7 
rho  =  1014 


The  elastic  constants  become 

rho  -  1014 
cll  =  2.11851*10^7 
cl2  =  9.69438*l(y6 
cl 3  =  1 29694*  m7 
c33  =  5.74394*10'^ 
c44  =  5.74575*10'^ 
c66  =  5.74539*10^ 

The  above  components  are  normalized  for  numerical  computations 

cll  =  3.6871  c44  =  3.7*c44 
cl2  =  1.68723  c44  =  1.7*c44 
cl3  =  225722  c44  =  23*c44 
c33  =  99.9686  c44  =  100*c44 
c66  =  0.999937  c44  =  c44 
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